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ABSTRACT 


I 

This  report  describes  two  computer  programs^  The  first  corrects 
the  orbital  elements  of  a  geodetic  satellite.  The  second  uses  residuals 
generated  by  the  first  program  to  correct  the  geocentric  coordinates  of  the 
observing  sensors  and  of  the  origin  of  the  datum  by  which  a  set  of  such 
sensors  are  connected. 


iv 


CONTENTS 


SECTION  PAGE 

1  INTRODUCTION  .  1 

2  ORBIT  CORRECTION  PROGRAM  .  A 

2.1  Theory .  5 

2.2  Formulation  .  7 

2.3  Flow  Charts  . .  32 

2.4  Input  and  Output  Formats .  49 

2.5  Operating  Procedure  .  58 

3  STATION  LOCATOR  PROGRAM 

3.1  Theory .  62 

3.2  Logical  Program .  64 

3.3  Flow  Diagrams .  76 

3.4  Input  and  Output  Formats .  98 

3.5  Operating  Procedure  .  103 

3.6  Experimentation .  104 


4 


MODIFICATIONS 


108 


AERONUTRONIC  U-1490 


INTRODUCTION 


SECTION  1 

INTRODUCTION 


Under  contract  AF  19(604) -7253 ,  Aeronutronic  has  developed  two 
computer  programs  to  perform  the  following  tasks; 

(1)  Orbit  Correction  Program 

To  compute  a  satellite  ephemeris  to  high  precision,  using 
optical,  radar  slant  range  and  doppler  observations.  High  precision  is 
understood  to  mean  that  errors  no  greater  than  20  feet  shall  arise  during 
any  one  day  portion  of  the  orbit  from  the  effects  of  zonal  harmonics  or 
computational  procedures  for  the  case  of  a  spherical  satellite  with  high 
mass  to  surface  area  ratio  in  a  circular  orbit  approximately  500  statute 
miles  above  the  earth's  surface.  This  error  does  not  include  errors  intro 
duced  by  air  drag,  radiation  pressure,  longitudinal  variations  of  gravity, 
observational  errors,  or  station  errors.  Air  drag,  however,  is  accounted 
for  by  a  model  atmosphere.  Radiation  pressure  and  longitudinal  variations 
of  gravity  can  be  introduced  by  modifications  of  the  program. 

(2)  Station  Locator  Program 

To  compute,  from  the  output  of  the  Orbit  Correction  Program 
a  geocentric  translation  vector  for  each  station  location  such  as  to  mini¬ 
mize  the  residuals  from  that  station  and,  for  each  datum,  a  geocentric 
datum  translation  vector,  which  minimizes  the  station  correction  vectors 
for  all  the  stations  on  that  datum.  This  program  includes  means  of  weight 
ing  observations  and  station  coordinates  according  to  criteria  determined 
externally  to  the  program. 
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The  decision  to  correct  the  station  location  in  a  separate 
least  squares  process  was  taken  partly  on  the  basis  of  expediency.  It 
was  thought  to  be  simpler  and  quicker  to  write  and  check  out  the  two 
programs.  Such  a  procedure  was  also  considered  to  provide  the  user  with 
greater  flexibility  in  the  use  of  the  programs.  The  following  practical 
advantages  may  be  cited  in  favor  of  separate  programs; 

(a)  The  matrices  are  kept  small  and  easily  manageable. 

(b)  The  orbit  computation  can  be  limited  to  a  few  revo¬ 
lutions  in  order  to  limit  the  effect  of  round  off 
errors.  Yet,  the  residuals  obtained  from  several 
such  orbital  arcs  can  be  used  to  correct  station 
location . 

(c)  The  observations  from  several  different  satellites 
could  be  combined  in  the  Station  Locator  Program. 

This  may  require  a  change  to  accomodate  several 
ephemer ides . 

(d)  A  different  form  of  orbit  computation  or  predic¬ 
tion  could  be  used  to  provide  input  to  the  Station 
Locator  Program. 

(e)  The  weighting  factors  may  be  changed  between  the 
programs  and  for  different  runs  of  the  Station 
Locator  Program.  This  advantage  vanishes  when  the 
evaluation  of  these  weights  is  formalized  and  can, 
therefore,  be  made  inside  the  program  or  programs. 


On  the  other  hand,  certain  theoretical  considerations  indicate 
advantages  of  simultaneous  correction  of  both  orbit  and  station  errors. 

The  principal  argument  is  that  what  is  desired  is  a  best  fit  to  the  obser¬ 
vations  of  the  total  physical  model.  This  model  includes  both  the  orbit 
and  the  observing  net.  Errors  in  station  coordinates  could  be  misinter¬ 
preted  by  the  Orbit  Correction  Program  as  errors  in  orbital  elements  and 
remaining  errors  in  orbital  elements  could  be  interpreted  as  station  errors 
by  the  Station  Locator  Program.  There  is  a  justifiable  fear  that  these 
misinterpretations  could  limit  the  accuracy  of  the  process.  That  is,  the 
programs  could  not  converge  to  as  close  an  approximation  to  the  true  model 
as  could  aprogram  which  corrects  all  quantities  simultaneously. 

The  same  problem  can  occur  in  the  combined  program  also.  That 
is,  the  confusion  of  error  source  can  be  attributed  to  a  similarity, 
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analytical  or  numerical,  of  the  partial  derivatives  which  express  the  re¬ 
lationship  between  the  error  and  its  possible  sources.  If  such  a  simi¬ 
larity  exists,  there  is  a  danger  that  the  large  matrix  may  be  nearly  sin¬ 
gular  and,  therefore,  unstable. 

None  of  the  above  argues  against  combining  the  two  programs  in 
such  a  way  that  separate  solutions  are  possible  and  that  much  of  the  inter¬ 
mediate  output  is  available  for  Inspection.  This  procedure  will  assure 
the  greatest  possible  efficiency  since  the  convergence  will  generally  be 
much  faster  in  the  combined  program. 

Other  desirable  modifications  will  be  discussed  in  Section  4. 
These  include  the  introduction  of  weighting  by  accuracy  of  observations 
in  the  Orbit  Correction  Program,  the  use  of  these  accuracies  to  give  un¬ 
certainties  in  the  elements,  the  inclusion  of  the  effects  of  additional 
perturbations  and  to  Improve  the  precision  of  the  model  for  atmospheric 
drag. 


This  report  is  Intended  principally  as  a  description  of  the  two 
programs  developed,  and  does  not  include  all  of  the  mathematical  develop¬ 
ment.  Those  who  wish  to  know  why  as  well  as  how  are  referred  for  theoret¬ 
ical  background  of  differential  correction  techniques  and  of  the  variation 
of  parameters  to  Aeronutronic  Publication  U-880,  Astrodvnamic  Analysis  for 
the  National  Space  Surveillance  Control  Center .  June  1960.  The  variation 
of  parameters  method  can  be  found  in  Vol.  1,  Appendix  3D,  "Efficient  Pre¬ 
cision  Orbit  Computation  Techniques,"  and  the  differential  correction  pro¬ 
cedure  is  shown  in  Vol.  II,  Appendix  4A,  "Differential  Correction." 
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SECTION  2 

ORBIT  CORRECTION  PROGRAM 


The  Orbit  Correction  Program  generates  satellite  epheiner ides  for 
the  computation  of  acquisition  coordinates  for  sensors  as  well  as  correct¬ 
ing  the  orbital  elements  of  the  satellite.  The  first  use  is  referred  to 
hereafter  as  simulation.  This  feature  is  incidental  to  the  tasks  of  this 
contract.  It  has  been  used  in  the  checkout  of  the  programs  and  may  prove 
valuable  also,  as  a  means  of  producing  look  angles,  scheduling  and  evalua¬ 
tion  of  sensor  sites. 

The  primary  purpose  of  the  program  is  to  correct  the  orbital 
elements  by  means  of  observations  of  slant  range,  range  rate,  right  ascen¬ 
sion  and  declination,  and/or  azimuth  and  elevation  angle.  After  the  last 
orbit  correction,  the  program  prepares  the  residuals  in  the  observations 
for  use  in  the  Station  Locator  Program  and  a  tape  containing  position  and 
velocity  vectors  for  the  satellite  at  the  time  of  each  observation. 

It  should  be  noted  that  the  observation  time  shall  be  the  time 
that  the  electromagentlc  radiation  left  the  satellite.  That  is,  the  time 
of  reception  of  the  signal  at  the  station  must  be  corrected  for  the  time 
it  takes  light  to  travel  from  the  satellite  to  the  station.  For  an  active 
radar  observation,  this  can  be  done  by  averaging  the  times  of  transmission 
and  reception  as  well  as  differencing  them,  but  with  passive  doppler  and 
angle  observations  one  requires  an  approximation  of  slant  range  in  order 
to  find  this  "light  time"  correction.  Neglect  of  this  effect  produces  a 
bias  of  about  25  parts  per  million  of  the  slant  range  in  the  direction  of 
satellite  motion.  This  is  a  bias  of  40  feet  at  a  slant  range  of  350  statute 
miles . 
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2 . 1  THEORY 


The  Orbit  Correction  Program  employs  two  computational  techniques 
which  speed  the  computation  without  loss  of  accuracy: 

(1)  The  satellite  motion  is  numerically  integrated 
by  the  variation  of  parameters  formulation, 
thereby  eliminating  from  the  integrands  the 
large  central  acceleration  term. 

(2)  The  differential  correction  uses  analytical  dif¬ 
ferential  expressions,  thereby  eliminating  the 
need  for  more  than  one  numerical  integration  over 
the  observation  period. 

By  developing  the  differential  equations  of  motion  in  terms  of 
parameters  which  remain  invariant  in  the  absence  of  perturbations  to  the 
two-body  motion,  the  dominant  two-body  term  is  supressed.  The  parameters 
employed  in  this  formulation  are: 

a  ■  eP  a  vector  defining  eccentricity  and 

perigee  location 

h  -/p  w  the  orbital  angular  momentum 

L  =  M  +  the  mean  longitude  of  the  object  at 
o  o  , 

some  epoch 

These  parameters  are  valid,  in  the  formulation  employed  here,  for  all  satel¬ 
lite  eccentricities  and  for  all  inclinations,  including  zero  in  either  case. 


The  differential  characteristics  of  a  slightly-perturbed  satellite 
orbit  are,  to  a  first  order.  Identical  to  those  of  the  osculating  orbit.  Thus 
the  cause-and -effect  linear  relationships  needed  for  differential  correction 
may  be  developed  analytically,  rather  than  by  the  alternate  "variant  calcula¬ 
tion"  procedure  where  neighboring  ephemerides  are  integrated,  in  each  case 
with  one  of  the  parameters  modified  by  a  small  amount.  The  parameters  em¬ 
ployed  in  the  development  of  the  differential  expressions  resemble  those  used 
in  the  variation-of-parameters  ephemeris  program,  i.e.: 


a  *  eP 


a 

U  -  M  +  <*) 
o  o 

n 

i 


a  vector  defining  eccentricity  and 
perigee  location 

semi-major  axis 

the  mean  argument  of  the  latitude  at 
some  epoch 

nodal  longitude 

Inclination 
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The  adoption  of  fl  and  1  to  describe  the  orbit  plane  orientation  and  of 

U  to  denote  initial  position,  restricts  the  development  to  non-equator ial 

orbits.  In  addition,  the  vector  a  is  described  in  terms  of  two  components 

in  the  orbit  plane  to  avoid  redundancy;  these  components  are  designated 

a  „  and  a  with  the  former  in  the  direction  of  the  ascending  node. 
xN  yN 

Any  observation  0^^  by  an  earth-fixed  observer  may  be  expressed 
in  terms  of  these  six  parameters  or  elements  describing  the  orbit  and  the 
time.  First  order  differential  expressions  relating  observation  and  para¬ 
meter  follow  from  the  leading  term  in  a  Taylor  expansion,  i.e., 

r 

where  the  X.  are  the  six  orbit  parameters.  Where  there  are  m  observations 
available  ti  define  the  six  parameters,  a  set  of  m  differential  expressions 
may  he  written;  in  matrix  form,  this  set  is 


(AO^)  -  (C^j)  (AX^) 


(2) 


where  is  the  m  x  6  matrix  t;lth  typical  element, 


<)  0, 


(3) 


The  (  A  Xj)  is  a  six  component  vector,  and  (^0^)  is  an  m  component  vector, 


The  elements  of  the  (C  ,)  matrix  are  the  partial  derivatives  of 
the  observed  quantities  with  respect  to  the  orbit  parameters  and  may  be 
determined  analytically  or  by  variant  trajectory  calculations,  wherein  the 
parameters  are  varied,  one  by  one,  and  the  resulting  changes  in  the  observa¬ 
tions  are  noted. 


If  there  are  more  observed  quantities  0,  than  parameters  X.,  that 
is,  when  m^6,  the  system  is  overdetermined  and  the  number  of  equations  may 
be  reduced  by  the  method  of  least  squares.  The  solution  takes  the  form 


(AXj)  -  (C^j)  ]  (C^j)^  (  AO^) 


(4) 
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where  the  -1  and  T  superscripts  denote  inverse  matrix  and  transpose  matrix, 
respectively.  The  bracketed  quantity  in  (4)  is  the  so-called  least  square 
matrix,  N: 


N  =  (C.j) 


(5) 


In  the  interest  of  achieving  efficiency, 

from  analytical  expressions,  which  are  detailed  in 
presenting  the  formulation  of  the  program. 


the  C..  are  evaluated 
ij 

the  following  section 


2 . 2  FORMULATION 

The  Orbit  Correction  Program  can  be  conveniently  divided  into 
three  parts,  one  concerned  with  the  correction  process,  another  concerned 
with  the  generation  of  different  types  of  ephemerides,  and,  most  important, 
that  portion  necessary  for  both,  (mainly  concerned  with  the  integration  of 
the  equations  of  motion).  These  parts  will  be  documented  below  in  the  in¬ 
verse  order. 


2.2.1  VARIATIONS  OF  THE  PARAMETERS 

The  major  portion  of  the  Orbit  Correction  Program  is  spent  in 
generating  position  and  velocity  of  the  satellite  at  successive  times  by 
numerically  integrating  the  equations  of  motion.  This  integration  process 
is  started  by  generating  the  initial  values  of  the  integrands  and  other 
quantities ; 


Given  a  „  ,  a  „  ,  h  ,  L  (subscript  o  indicates  epoch) ,  the  fol- 
xN  yN  ^  o 
o  o 


lowing  procedure  is  common  to  both  the  simulation  and  differential  correc¬ 
tion  portions  of  the  program: 


(a) 


Compute  Greenwich  sidereal  time  at  epoch. 
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(in  degrees)  =  D  (©-360  )  +  ©f  +© 
gr  '  ®  gr 


where  D  is  epoch  day  number,  f  is  the  fraction  of 
a  day  elapsed  from  start  of  epoch  day  to  epoch, 
0^  is  Greenwich  s.idereal  time  at  the  beginning 
of  §poch  year,  and  G  ■=  360?9856472,  the  rotation 
rate  of  the  earth  in  degrees  per  mean  solar  day. 


(b)  Compute  the  semi-latus  rectum: 


p  “  h  ■  h 
'^0—0—0 


(c)  Compute  the  orientation  vector,  W^i 
W  = 

/•/  P, 


(d)  Compute  the  orientation  angles,  i^ ,  .  >^,and  U^: 


sin  i 


cos  i  =  W 
o  z 


,  sin  X 
-1  /  o  , 

tan  ( - - — ) 

cos  i 


0  ^  i  d  TT 
o 


W 


>in  SI 


°/ 


sin  i 


-W 


cos 


a 


7-  ■ 

/sin  1 


n 


,  sin  fl 
►  "1  /■  o 

tan  ( 


cos 


0  <  z-rr 
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U  =  L  +.. i  If  W  <0  (retrograde  motion) 
0  0  0  2 

o 


U=L-i^ifW  >0  (direct  motion) 
o  o  o  z 

o 

(e)  Compute  the  equatorial  coordinates  of  a^ : 

a  =  a  „  cos  -  cos  i  a  „  sin 

X  xN  o  o  vN  o 

o  o  o 

a  /  a  =a„  sinj^  +  cos  i  a  „  cos  /i 

— o  \  y  xN  o  o  yN  o 

'  o  o  o 


a  =  a  „  sin  1 
2  yN  o 

o  o 


(f)  Compute  the  eccentricity,  e  ,  and  the  semi-major 

axis,  a  : 
o 


2  2  2  ^  2 

e  «  a  +  a  +  a 
o  X  y  2 

o  'o  0 


°/i-< 


•o  ■  'Vl-e  ^ 

o 


Also  common  to  both  simulation  and  differential  correction  is  the 
numerical  integration.  The  equations  to  be  integrated  are  of  the  form: 


dyf 

{jt  1^2’  ^6’  ^7^  i  “  1 »  2 . 7 


where  the  y.  equal  a,a,a,h,h,h,  and  L. 
i  X  y  z  X  y  z 

The  numerical  integration  scheme  used  is  based  on  the  following  fourth  order 
Runge-Kutta  method: 


n+1  n  .  At 
^i  — 


(K,  +  2K.,  + 


2K_  +  K, 


) 
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where  •«  (t  •  Vy) 


2 

K„ 


K  =  £  (t*'  +  — ^  -  y”  +  — ^ 

3.  i  '■  2  ’  2 

L 


n  ^27  . 
...y;  +  -^) 


=  f.  (t"  +  At.  y'l  +  K3^,  ....  y^  + 


da  dh 

As  can  be  seen,  it  is  necessary  to  compute  — 7~~y  - 1  and  —f 

dt  dt  d 

from  a,  h,  and  L,  several  times  at  each  integration  step.  The  first  step 
is  to  compute  position  and  velocity,  r,  r  from  a,  h,  and  L. 


Given  a_,  h,  and  L  at  some  time  t,  the  following  procedure  is  used 
to  derive  position,  r,  and  velocity,  jr.  Note  that  whenever  the  anomalies 
V  ,  E,  and  M  are  used,  they  appear  either  in  a  sum  with  CO  or  in  products 
with  the  coefficient  e.  Thus,  no  indeterminacy  exists  for  zero  eccentricity. 
In  its  present  form,  precisely  equatorial  orbits  cannot  be  integrated,  since 
the  ascending  node  is  employed  as  a  reference  direction. 


(a) 


Compute  p,  e,  a,  n: 

2  2  2 

p=h’h  =  h  +h  +h 
r  _  _  X  y  z 


2  2  2  2 
e  =a'a=a  +  a  +a 
—  —  X  y  z 


n 


3/2 


where  =  .07436574 


(b) 


Compute  the  orientation  vectors  W,  M,  and  N: 


(note  W  =  cos  i) 
z 


M  =  /  1  -  =  sin  i 

z  V  z 


N  = 

X 


cos 
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M  ■=  N  W  “  cos  n  cos  1 
y  X  2 


>ln  S2 


M  =  '  N  W  =  -  sinJ7  cos  i 
X  y  z 


N  =0 
z 


(c)  Compute  the  components  of  a  in  the  orbit  plane,  a^^j  and  a^^ 


xN  =  a  ■  N 


a  =  a  •  M 

yN  —  — 


(d)  Compute  the  orientation  angles ,  and  U: 

Cl  =  tan  ^  ^  ^  O  ^  Cl  <  2'n' 

U=L+X2  ifW^<  0  (retrograde  motion) 

U=L-«riifW  >-0  (direct  motion) 

z 

(e)  Solve  Kepler's  equation  for  E  +  CO  by  iteration,  using 
a  first  guess  of  U  (mod  2Tf  ): 

E+CO=U+a„  sin  (E  +  CO)  -  a  „  cos  (E  +  CO  ) 
xN  yN 

(f)  Compute  It  and  £; 

e  cos  E  =  a^jj  cos  (E  +cO)  +  sin  (E  +60) 
e  sin  E  =  a^  sin  (E  +60)  -  a^^^  cos  (E  +6<J) 


r  =  a  (1  -  e  cos  E) 


e  sin  E 
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rv  /l  -  e^' 

COS  u  =  "f"!"  (E+dO  )  - 

sin  (E+6>J)  - 

U  =  cos u  N  +  sinu  M 
V  «=  -  sin  II  N  +  cos  u  M 

r  =  rU 

£  •=  fU  +  rvV 

From  the  position  and  velocity,  it  is  possible  to  compute  the  per¬ 
turbative  accelerations,  specifically  the  bulge  perturbation,  f„: 

j'(5U^  -  1)  +  ^  h'(7U^  -  3) 
z  /  z 

r 

/  ?  L  S  7  L 

■  K  (42U  -  63U  -3)  +  — 5-^  J,  a  u.  (5-30  U  +  33U  ) 

z  ^  g^8  5  er 

.  I 


'B  5 
r 


6r 


•  '  _  y/ 


sin  u 
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4  =  -f  -  3  )  ^ 

r  5r 


+  k'^(-  63U^  +  70  -  15)  + 

^  •  Z  Z  O' 

6r  8r 


4  2 

U  -  lOU  +  1) 
z  z 

J,  a^/«  (-5  +  105U^ 

5  e  z 


-  315  +  231  U^) 

z  z 


where  U  =  -  and  j'"  =  .001  623  41 

z  r 

h'  =  -.000  006 

k''  =  .000  009  09 

J,  /l=  -.000  000  2 

5  e/ 

At  each  point  during  the  integration  at  which  the  derivatives 
dL/dt,  dh/dt ,  d^/dt,  are  evaluated,  the  perturbative  accelerations  due 
to  drag  y^,  z^,  must  be  evaluated  and  added  to  the  bulge  accelera¬ 
tions,  x',  y',  z',  to  obtain  the  total  perturbative  accelerations,  x'  y' 
,  .\  b  b  b 
and  z. 

Given  ^  and  and  tabulated  values  of  the  density  ratio  and 
atmospheric  molecular  weight  versus  altitude: 

(1)  Compute  the  relative  air  speed  vector,: 


=  X  +  y  e 
Vy  =  y  -  X  0 

y  =  z 

'z 
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Also  compute  the  magnitude  of  the  relative 
air  speed  vector 


V 


6  In  the  above  equations.  Is  the  angular  rota¬ 
tion  rate  of  the  Earth. 

(2)  Compute  the  altitude  above  the  oblate  spheroid 
in  Earth  radii: 


H  =  r  -  1  -  I  f^  +  (f  +  I  (f)^, 


where  f  la  the  flattening  of  the  Earth  *  1/298.3. 

(3)  Look  up  the  density  ratio  0”=  CfCH)  and  the 
molecular  weight  M  =  M  (H)  from  the  tabulated 
data  (see  Appendix  C  for  the  programmed  tables), 
and  calculate  the  atmospheric  density  : 


o 


where  P  is  the  sea  level  value  of  the  atmos¬ 
pheric  density  =  0.001225  gm/cm^. 
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(4)  Compute  the  skin  temperature  of  the  vehicle; 


pc  (V  , 

S 


1/4 


where 


(f'  =  Stefan-Boltzmann  constant  =  5.672  x  10  , 


and 

C  =  emisslvity  of  the  satellite  =  .9. 

(5)  Compute  the  drag  coefficient,  C^,  by  first 
computing  the  auxiliary  quantity 

.  6.972  X  lO^Vd 

C  -  t  ■  »  > 

C„  VQA  T  ) 

Do  y  e  s 


and  then 


.  -ctr 

C„  =  C„  (1  +  1.1739130  e  )  , 
D  Do 


Cj.  being  a  reference  value  of  the  drag  coef¬ 
ficient  =  0.92. 
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(6)  Compute  the  drag  terms: 


A 


X— ►y,  z, 


where  d  is  the  diameter  of  the  satellite, 
m  Is  Its  mass,  and  K  Is  a  constant  rela¬ 
ting  the  units.  The  program  value  of 
K  =  2.504742  x  10^. 


Finally  we  are  ready  to  calculate  the  Inte¬ 
grands  da^/dt,  dh/dt ,  and  dL/dt: 

(a)  Determine  the  total  perturbations 


.\ 

r  = 


% 


(b)  Compute: 

r  =  r  • 

.  .\  . 

s  s  =  £  ■  £ 

r  r  =  r  •  r 

and  the  auxiliary  quantities: 
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D  - 


I 


•\ 

D  = 


(c) 


Compute  n,  rl^, 


\ 

n 


>  A 

rb  =  W  •  f 


t-  _ ^ 

(1  + 


^  -  D£  - 


e2[  =  W  X  a 


-  e  v=  e2  •  ^ 


and  L: 
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\ 

L  - 


2D’  .  eV 

yr  1  +  "/l  - 


(d)  Compute  h: 

.\ 

\  £  X  I. 

-'W 


(e)  Compute  the  derivatives: 


dt 


k  l"  +  n 
e 


d£ 

dt 


k  a' 
e” 


dh 

dt 


k  h' 


2.2.2  EPHEMERIDES 

The  word  ephemerls  is  defined  as  a  table  of  the  positions  of 
celestial  bodies  at  regular  Intervals  of  time.  As  used  here  In  refer¬ 
ence  to  the  Earth  satellite  ephemerldes,  the  Intervals  of  time  are 
Integration  steps.  The  ephemerldes  are  of  three  kinds. 
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The  first  type  of  ephemeris  shows  geocentric  position  and 
velocity  vectors.  No  additional  fomulation  is  needed  to  produce 
this . 


The  second  type  of  ephemeris  shows  the  position  of  the 
satellite  as  latitude,  longitude,  and  height  above  the  geoid.  This 
is  the  subsatellite  track. 

The  third  type  of  ephemeris  shows  acquisition  coordinates 
for  a  number  of  designated  sensors.  These  "look  angles"  show  the 
coordinates  at  which  the  satellite  would  appear  to  the  station  if  it 
had  ::tarted  with  the  given  boundary  conditions  at  the  epoch.  These 
acquisition  coordinates  are  produced  for  every  integration  step  at 
which  the  satellite  is  above  the  horizon  of  the  designated  station. 


a .  Computation  of  Sub-Satellite  Track; 


If  the  option  to  compute  a  sub-satellite  track  consisting 
of  latitude,  <f>  ,  East  longitude,  Ag,and  height  above  the  earth,  H, 
is  chosen,  then  these  quantities  are  computed  at  each  point  of  the 
ephemeris  as  follows: 


<P 


tan 


-1 


U 


(1-f)' 


90°^  ^  <  90*^ 


where  f  =  ■  —  is  the  flattening  of  the  Earth 

270  «  J 


A  (in  degrees)  =  ©  -  ©  (t-t  )  -  ©  A-,  ^  0° 

E  o  gr^  E 
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where  0  =  tan 


-1 


0°  f  e  <  360° 


and  0  =  .250  684  48 


is  the  rotation  rate  of  the  Earth  In  degrees 
per  solar  minute. 


H  (in  Earth  radii)  =  r  -  1  +  (^  +  f)  | 

i  z  i.  z 


b .  Simulation  of  Acquisition  Coordinates; 

This  part  of  the  program  simulates  "observations"  of  the  satel¬ 
lite  in  the  orbit  specified  by  the  input  parameters.  Given  0,  the  lati¬ 
tude  of  a  station  in  degrees,  the  East  longitude  of  the  station  in 

degrees,  H,  the  height  of  the  station  above  sea  level  in  meters,  and  0  , 

Greenwich  Sidereal  Time  at  the  initial  time,  t  ,  in  degrees,  the  follow¬ 
ing  procedure  computes  0(,  ,  A,  h,  ^  ,  and ?or  time  t: 

(1)  Convert  0,  X  ,  and  0  to  radians  and  H 
to  Earth  radii 


(3) 

(4) 


Compute: 

C  =  (1  -  sin^  0)’ 

^  ^  ^  298.3  ° 

of  the  Earth. 

S  =  C  (1  -  €^) 

Compute: 

0  =  .004  375  269  1  (t-t  )  +  0 

o  gr 


Compute  the  station  vector,  R: 


the  flattening 
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X  ■  -  (C  +  H)  cos  (j)  cos  0 

Y  ■  -  (C  +  H)  cos  p  sin  © 

Z  -  -  (S  +  H)  sin  ^ 


(5) 


Compute  the  slant  range 


^  +  R 

P  “  ^ p  ‘ JQ  -^(x  +X)^  +  (y  +Y)^  +  (z  +  Zp  ' 

(6)  Compute  azimuth,  A,  and  elevation  angle,  h,  of 
object  as  "seen"  from  the  observation  station: 


L  .  (x+X)  cos©sinj^  -t-  (v4-Y)  sin©sin^  -  (z+Z)  cosp 

xh  yO 


T  -  (x+X)  sinQ  -f  (v+Y)  cos© 

yh  "  p 


'zh 


(x+X)  cos  ©  cos  ^  -f  (yfY)  sin©  cos^  +  (z+Z)  sin^ 


tan 


P 

,  -  TT/2  ^  h  ^T/2 


If  h  ■<  0  (i.e,,  the  object  is  below  the  local  horizon) 
the  output  for  this  time  is  omitted. 


A  •>  tan 


^  A  <  27r 


The  quadrant  is  determined  from  an  examination  of  the 
sign  of  the  numerator  and  denominator. 

(7)  Compute  topocentric  right  ascension,  Ql  •  declina¬ 
tion,  5  • 
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(L  ,  L  ,  L  ) 
X  y  z 


-ir/2^S^'iT'/2 


a  = 


0  ^  a  <  2Tr 


The  quadrant  is  determined  from  sign  of  numerator 
and  denominator . 


(8)  Compute  the  slant  range  rate,^  : 
X  =  -  Y© 

Y  -  X©  )  ©  -  0.058  834  47 


yO  "  i  +  i  (X,  Y,  Z) 

y©  =  L  (x  +  X)  +  Ly  (y  +  Y)  +  (z  +  Z) 


2.2.3  DIFFERENTIAL  CORRECTION  OF  ORBITAL  ELEMENTS: 

This  part  of  the  program  relates  residuals  in  the  observations  at 

time,  t,  to  corrections  to  be  applied  to  the  initial  orbital  parameters  at 

time  t  . 
o 


The  procedure  calculates  the  orbital  parameters  and  quantities 
associated  with  the  station  coordinates  at  the  observation  time.  It  com¬ 
bines  these  quantities  to  obtain  the  coefficients  of  the  linear  relation¬ 
ships  relating  residuals  to  any  combination  of  Aa  /  a  ,  Aa  „  ,  Aa  „  , 

o  o  xN  yN 

represents  the  observations  to  determine  the  res?duals 
Finally,  the  corrections  to  be  applied  to  the  parameters  are  determined  by 
solving  the  (usually  overdetermined)  system  of  linear  correction  equations. 
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of  the  observing 
station 


a .  Forming  the  Linear  Correction  Equations 
Given: 

(1)  ,  the  latitude  in  degrees 

(2)  X  ,  the  longitude  in  degrees 

£ 

(3)  H,  meters  above  sea  level 

Greenwich  Sidereal  time  at  epoch  in  degrees 

(5)  t,  the  time  of  the  observation  in  minutes,  and 

(6)  one  or  more  observed  quantities  at  this  time, 

the  following  will  compute  one  line  of  the  above-mentioned  system  for  each 
observed  quantity. 

(1)  Repeat  steps  1  through  4  under  simulation 

to  obtain  &  and  ^  in  radians  and  the  station 
vector  R. 

(2)  Compute  the  coefficients  R  and  U  where: 

R^  =  (a^/r  ^  e  sin  E 


‘^a  “  ^  'J  -V  \ 

‘  [“xN  "  +CO)] 

V  “  [“yN  ■ 


^  f  (U  -  U^) 
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a 

=(k,  • 

U)  R 
—  a 

•  V) 

u 

a 

xN 

■  £>  'xN 

•  V) 

^xN 

C  A 

■  2>  > 

•  V) 

^yN 

"au 

=  (L, 

■  U)  R 
—  u 

+  (i^ 

■  V ) 

U 

u 

o 
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■=  ■  V)  r  cob  i  -  (1^  •  W)  r  sin  1 


C  y.  .  =  (L  •  W)  r  s  in  u 
Ai  -c  - 


Enter  the  following  linear  correction  equa¬ 
tion  into  the  system  of  such  equations: 


A.  -f  ^  'a.  ^ 

I  —  o  xN  '■ 

a 


yN 


(5)  If  A,  azimuth,  and  h,  elevation  angle  are 
observed,  compute; 


=  sin  ^  cos  © 
Sy  =  sin  JZS  sin  & 


cos 


sin  0 
cos  Q 
0 


=  cos  p  cos  © 


cos  u 
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L  .  ■=  -  cos  A  cos  h 
xh 


"yh 


zh 


sin  A  cos  h 
sin  h 


h 


xh 


yh 


zh 


sin  A 

cos  A 
0 


D  .  =  ccos  A  sin  h 
xh 


D  ,  =  — sin  A  sin  h  V  D. 
yh  /  -h 


D  ,  “  cos  h 
zh 


L.*L,S+L.E+L.Z 
-obs  xh  —  yh  —  zh  — 

A,*A,S+A.E+A,Z 
-obs  xh  —  yh  —  zh— 

'd.-d.s+d.e  +  d.z 

-obs  xh  —  yh  —  zh  — 


Compute  AL  ■» 


Form  the  coefficients  as  in  (4)  with  A  . 
replacing  L  and  enter  the  following  linear 
correction  equatiion  into  the  system  of  such 
equations: 
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(6) 


/^c^bs  •  ^ 


a  a^ 


^  Aa 

xN  o 


+  C.  Aa„+C..  „Au  +  C  . ,  Ai'i  +C..Ai 
Aa  „  vN  AU  o  A‘,'  o  Ai  o 
yN  o  o 


Again  form  the  coefficients  as  in  (4) ,  this  time 
with  D^jjg  replacing  L  ,  and  enter  the  following 
linear  correction  equation  into  the  system  of 
such  equations: 


P 


D  ,  • 
c-obs 


Al  =  c 


Aa  a 

-  o 

a 


+  C. 


xN 


xN 


+  C 


yN 


yN, 


+  c 


AU 


A  u. 


+  C 


A  i 


If  CH.  ,  topocentric  right  ascension,  and  O  ,  topo- 
centric  declination  are  observed,  compute: 


^bs 


A^  =  -sin  Oi 


A  =  cos  (X 

y 


A  K 
^bs 


A  =  0 

z 
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D  =  -  sin  o  cos  CX 

X 


D  “  -  sin  0  sinO- 

y 


D  . 
-obs 


D  =  cos  6 
z 


Compute  AL  *  ^iobs”  ^ 


Form  the  coefficients  and  compute  the  linear 

correction  equations  as  in  ^)  ,  substituting 

A  .  for  'K  ,  and  D  .  for  D  , 

—obs  —obs  —obs  —obs 

(7)  If ^  ,  the  slant  range  rate,  is  observed, 
compute; 


X  =  -  Y0 

Y  =  xe 

Z  =  0 


R  where  B  =  0.058  834  47 


A 

.4 


=  r  +  R 


A 


ex^  ■  a  (e  cos  E  -  e  ) 


/ 


eyu>  «aJl-e  esinE 


rv 

r 


Compute  the  coefficients  R  and  U  where: 


3/2  ,  3 

a  exw/r 


r/2  -  3/2  (U  -  U  )  R 
o  u 
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^xN  “  ^  sin  (E  +CJ)-  e  sin  E  - 

V  "  cos  (E  +6J)-  a^j^  e  sin  E  + 

IT  r~r  ^ 

-  -/yt/t  a  ey^  /r 

-T  -  I  -  "o>  "u 

“x»  ■  /  1  -  ■  l^cos  (E  +c^)  -  +  ^)j 

(8)  Form  the  coefficients: 

■  (t,  ■  !D[^,  («,  -  V  +  (^  ■  E)  E. 

a 

+  (1,  •  'D\f>,  «i.  -  i  ".)-/Oc>'a>  <;4  ■  »  "a 

<\n-  ”  “aN’-yOc^N]  *Yl  ■  a>  "xE 
+  <Ea  •  S)1)0,  i  "xl.>-/Oc%]  *  % 

<ic  •  2)  |/>e  <V-  -  >)-  ^Eys]  +(/2c  ■  E) 

■  2)[/o,  (V  r  v>-  A  v]  *  A '  -’  V 

C  0^  ■  <!x  ■  2)[/»c  A  -  E  “u>-Ac  'u]  ^  <A  ■  -'  "u 

+  (E,  •  V)  A  f  I  ii„)-  u;^+  •  V) 
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■  -y^c  i  +  (kc  ■  ‘=08 

+  (yOc  ■  ‘=08  i  •  W)  sin  i^  (rv  sin  u-v  cos  u) 

"A  r  cos  o  J  -  •  W)  r  sin  i  cos  u 

/°c 


=  (I^  •  M)|yC).  ‘=08  o  +  f  sin  u)  -  /0_  r  sin  u  j 


J 


Compu 


+  (  O  ■  W)  r  sin  u 
^c  - 


Enter  the  following  linear  correction  equation  into 
the  system  of  such  equations; 


A®. 


P 


/ 


+c  .  Aa.„+C^„  A  ^..M  ^AT?  Au 


^AyO  ^  a^  '  ”  Aa  *“*“xN  '  “Aa  „ ‘~'“yN  '  'AU  o 


xN 


yN 


♦■tva  <1/1 


b .  Computing  the  Corrected  ,  a^j^  ,  a  „  ,  h^  ,  h  ,  h^ 

0^0  O  o 


Let  ZL  C..  AX,  =  A  0,  i=l,2,  3,  ...  represent  all  of  the 
j  =  1  =-J  J  ^ 

linear  correction  equations  (i.e  ,  the  C.,,  are  the  coefficients,  theAX.'s 

ij's  j 

are  the  corrections  to  the  orbital  parameters  at  time  t  ,  the  AO.  being  cor¬ 
rected).  The  following  matrix  equation  is  solved  to  give  the  corrections, 
in  a  least  squares  sense,  to  the  orbital  parameters  at  time  t^. 
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^12 


-<^11  ^12 


•E"i2  "iN  ^^2  Z"i2 


'^12  ^iN 


AX^  Zc^N  Ao. 


These  corrections  are  applied  as  follows  (a  prime  means  that  the 
element  is  a  corrected  element): 


®  XT  ^  ®  M 

xN  xN  xN 

o  o  o 


a  \t  “  *  XT  ^  3  XT 

y^o  y^o  y^o 


a*  =■  a  (1  + -  ) 

o  o  a 

o 


i'  -  i  +  At  0  £  i  <11 

O  O  O  O'' 


=  L  +  AU  -  if  cos  i^  0 

^  ^  ^  Q  Q 


o  o  o 


l'  -  l  +  Au  +  Afl  if  cos  i*  >  0 
o  o  o  o  o 


W  ^  =»  sin  i  ^  sin 


W  ^  -  sin  1^  cos 

y  o  o 


=  cos  i* 
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I  2  I  2  I  2 

xN  ®yN 
o  ^  o 


p'  =  a  (1  -  e' 
o  o 


h 

-t) 


'  -  Jpt^  W 


/ 
— o 


a' 


/  .  /  .  r\/  / 

- .  o'k  “  a  » 


-  a  „  cos  1  sin 
vN  o 

^  o 


a  cos  cos 
o 


cos  /i' 
xN  o 

o 


yN 


/  / 

a  „  sin  i 
yN  o 

■'  o 


A'  +a' 


„  sinoA 
xN  o 

o 


2.3  FLOW  CHARTS 


The  following  diagrams  indicate  the  computational  procedures  in 
diagramatic  form.  The  first  two  pages  illustrate  the  general  flow  of  compu¬ 
tet  ton.  The  succeeding  pages  contain  a  more  detailed  view  of  the  program. 


In  the  diagrams,  ovals  indicate  subroutines  of  the  program.  They 
are  identified  by  their  symbolic  location  in  the  oval.  The  rectangle  blocks 
contain  other  computations,  bookkeeping  instructions,  etc.  Circles  indicate 
connectors.  In  the  detailed  flow  the  circles  usually  show  the  symbolic  loca¬ 
tion  names.  The  single  letters  do  not  correspond  to  actual  program  locations 
Diamonds  are  used  for  logical  decisions  made  in  the  program.  The  conditions 
causing  the  program  to  proceed  In  the  different  directions  are  given  beside 
the  appropriate  arrow. 
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f  E-*vter> 


various  vo^,-w 
■'•V«  V-Kte^rtfcvow. 

^  of 

ftOcK 

<l‘v^ev(!tV>f  V9  -wo^'te 
tV<  CM'ntL  wiilWfi'nt 


to) 


rfr  re^vk\aor 
\>\^ttr«^'«on.  pt. 

Wt  lotypvii  tW 
t  i-wR  Wfc  ^p  M;  li 
conrMujEe  restctft  el+t.  • 
'Tvert.  W  rettw^  «.ltii. 
'4Sr«  unri«uXt/i  VttX 
tv*mt  re<Tort  at  Oud 
VMtt  rt&<u^ 


YVocefts  oWsfirvoticms, 
wri'ti'n^  04  W\omi  oh 
P  W  o.  p  ols.j 
w  ««"<,>  oU.^ 
fer  oo\  oW, 
^cu*A  ^  Wo.  p  ol«., 
p\v*4  t-"t,oc*A  t  S»1' 


oi-^d  taps 
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OvjuVpULt  *-esi(iuLa\s 

last  oIa-  <»^to  ta^e  s 

■9or  csW-lvne  ^a»cWii«4 
(x-vui  prl-Kfe  yw^ .  ° 


Write  "titjC  oA  -ta^ej 
t  tVvkS  is  durvwA 
•vUe  iMkSs  "tUru. 
obs.  i-»il6^in^ 
•t.V\i>.  eaV<e>^ris  "* 

I  Co'wv.jxitfi  Rok. ,  R*»ij 
^vjw,Ra,UQ^  > 

■i^k-i^-a]|^  I  Jt  i>lft&- 

p>iter  CoEFF  suWrBut  ime 
tb  covtfvita  otui  n»tt)  ^ 

coe(fs.  W  Le»^ 

St^Okores  yiotrix. 


ROOk^  '^iX  tabs, 
•tn>->n  tope 


TteokA  e^W-wkeris  ■to^e, 
C3^  So  'tWkli'  oV)S.  tzvwve. 
\S  VxMcyvded  Wvj  3  ti-NUSs 
o>v,  eitUer  %{(!€• 


Output  tWi  residaoi 
■W  oW-liifie  pri'nti'VkA  I 
uii-Hv  QSteri^  cct  *' 
tVke  ri^Vvt. 


P>Ai  -iC^uAre  cst 
reSktduai  Ukvltj)  rwisj 
SVXM.V^f^  vO  A 


ESKter  LCokSt  SSt.vjkA.res 
Routwive  ^  UtSrMvUiibe. 
rt)vO  CoePrs.  vu,ti 
Leost  '2<ykS>.ine5  fAt^vix 


Do  C.  ^ 


i-ntcrpoVoiticTTv, 

at  oU.  tVvAe 


v-ntCr^i 

L-,t^jVs 


PriVvt  ot4-Vv>v^ 
■nevib  G,Vts  a«vA 


EyXir  KokSt  Ssj-oans  . 

rovJk.t:i-»ve.  Ti  So\ue,  j 
UeokSV  SaoBfes  Hi{rn<  J  A  \ 
HvTccrrs.  Ti 

bWe  o-it^ _  <woftr»». 

CovnpvXtfe  r-wis  of 
rccn^  Oa\A  ocn^t 

iresiidvikAW  0k«^ 

W  TttM^  wife  resiiiiuj. 

fViYvt  oft'Vvwe  -vWtq 
rvvAts  <>\v)kS  tUe 
cuarrs 


OvitpyCL  restart 
U/CTor-v^OCtiov  oVklo 
tok.^  s  ^-or  ott- 
l»se  prt-rvti-K^ 


‘SdT  up  residvxcU 
rejectvorn  eriterux. 
^or  vveact  paO-SS 
tVyru.  ttvE  oWa. 


Vslrokpu^  "tape.  S 
■tor  otr-  We 
pa^TcWrn^  cxxvd 
pri-rtv-wA  ■ 


Write  oTfi  -Vletovori^d 
'ftLLS  VOeuLTMAT 
ttoDS*' 


E>vter  KVtSe  Sata- 
rokJkt  \'*ve  IS  ejowvpotfi 
Y:,t  ^as  Vnteirvwaiiok^el 
(^ucoiaities  ctt.  dU.tiWe 


WA  ccrre^vcnrvi  Xi  ^*t, 
TjevM  (Xo.OHav^ajjM,, Oo, 
Too  >  ojrA  io  • 

Eyiter  it'fltJs©  Siulorwjiti'yvG 
to  Cap-wftpute.  'Vkeio  Co.C, 


VOrVte  Grror'rwaVe 
earn  ftexouariter 
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|Wr;tt 
jtcLpes  3,  (.3 

VO 


WCT  i .  «TAeri 
/^uii-nd  \a^  o  \  ll&ui  > 

V«  t 

^cardB^  ; 
\^'’oca9S  (xmA  «tare 


BEGIW 


(.X'AtHSG^ 


VSL  *■'8^ 
r^T-* - 


CEFUhCy 


/^PyNCH.COS  > 
[  '3«ec.i<ki  X  oA  Aoitx. 
yw'tc.'t  c)w.  W- 
NflK-U*«_auwncWa4/ 


vx)  — *  SKtPSAV 
*TJK|  «J<OlT 
TMP  n>»K*->»3 


I  noTART  I 

\«o/w»TnicrtN. 

r-3^  V  — »  RSTMTC? 

H  ^<0 

—^CT^cq,  4 

VOXA^  — p  KXXX 
0  — VKTXPfcfc 
0  — ►  MT«f  fc»  j 

p/o  — »Xl4Tv«)0 

^nv  C>nVlD»C 


0*.»T1&« 


'iesy/^ Lt-MTMoNw 
CUTWD 

I - * 

.— *C4>moo+<^ 

'■V  -tV^re  *»  (X  ., 
rwtojrt  tw»\e  ,\i 
-nbt  f /o-*oJTwi»^ 

V  , 

T^Xuix-nA  to^  ^ 
r/o  — ♦  TPr& 

O  — V  KTAP%C 

— ♦GRTfcV 
(.^avvnai  -•Gilirsa 
v.o  — ►  0»*CWT 
F/b  — •  oiACTa. 
*TX»AS  L^TCruOt  L0N» 
nooE  altitup^^ 
CWKI  WrwML  CMSO-t-M- 


Cctmif^A 

oA  cfocVv* 

T 

V^— ^tTP 
U  — ►vJX,v4Y,Wl 

»W\  jL  - *StT>*t 

Ri'A  A  — *  SXKiO 
A  — •  CxsSo 
Co*u.  — *Cost 
o«x<Tr  —sYtAtv 
oiA<xTr  — > 

VC  Mo  DC  xV 
Se  — *-  Ax,  AY,  KX 
x-t*-  — >•  cs<a 
a.  — >  R  t,  AXV 


V  - •NCMTXV 

SS  — *  V.XHS 

r/o  — »coo^^T 

V  — •tPCcNT 

* 

l-XM  - *R31N1 

Ui*  - vRYRV 

TO^  - »«.«5Tt) 

A  — *RKsrn>+3 

>«  - vU.KVTO'^V 

X<  — >RK3Tt>t-S 

L - »RK3T0  Vt 

- >t.V.%T0^’1 

4  — VRK3T04-? 

tW*u  Kkstq  ■»■  lo 


0  — *  HTRPta  I 

,  i 

*J>\P  OERtV 
JMP  CATRv"— 
‘RXCNT  TIB,  'rTy», 
_ PTVY*^-^Q. 

/^hge.aomgeN 

[  x^taAmte.  ) 
V  e^Vvewiens  y 

'^1  ^ 
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CONTKOI.  ROUTINE  FOR 


CPHCMCRIS  INTE6RAT 


I 


(cw 


f*  -•EPl\tM'j»CVITIj3 


36 


|(t«-no-*-v  iVu. 

S.KC,TovJ.6',  = 

— ► 

Ttsic*snu. 
Hvru.  ePHtM 
•vLMCviTR.u'i+C 


I 


>Q.  CNThrtiwt 
-»NTftPfcC 

TuTA  0V\  iLVO 

— >xnkT 


LkJC.MTTRl'^  +  1 
— 5  NCHTRU 

r 


tCOU»4T^-*-V 

- ?><jOOt4T 

CfltMTROi-H. 

— *MC>JTtL. 
’Tvuv  St.^.0 


Te-st 


mmm 

( 

HtPlO  \ 
PrlwX  90^) 

5  — » \.X»46 
Cx?&CMT)+V 

— »XPG.CNT 

T— 

L. 

(^UX'te'i  -V  V 

- ►  l-XMfc 

\0ri%e  lEPHtM 
tV«ru. 

0V'^:A^  3 

tKfTBP63)-»l. 
_»NTt\9t3 

UiJ  wrw  o-A  V»i.V 
t— »  mC’^TR'- 


PRXNfT 
Xe  >  ^  er»\.”'iU^ 


C.X  -tV'Jk.  xvsp 
CMT-jjOtal  -tKra 

c-vjTvuo^at 


i  ^ 

^p^-vuiKt  No 

sjt.  we  be*.0Tvi>^^ 

iq 

^  it  «  Lte-t,> 

1 

-Lx -Vi  -  fce-t 

— ' 

^0  -»ciJruJD+vi 

it 

<^An-ujo»jo 

tV'x  tjjr<jc>  +Hijj 

CYV.  tape  (d 


^  ■ 

^OBit '■cs^Tfittd 

V.lr«AaF)  — • 

t*JTvjOV>-»1.3 

h 

tvxfv.  i_OToe+l'» 

l\JOO  — aCJJTViX)  VVM 


lit 
Ri. 

•Set  fi.t^ak'C 
^  ’Ofett  »««W 
JFttr 

rfcsnvrTtiMt 
•»0  -iKjJTjtm.* 
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(NTAre«)-*ISAVt6 
60  -*>  UINC 
I  —  tricuT 
"  STA  TIMC-MIN 
RANAE-IUBIt  RHOPOT 
-  KAO/kCHIM  RT  ASCCN 
_t>e6  ZlECi.-Pe6 
AZIMUTH  -PC6 
ALTITUDE- PE6" 

— »CHCP 
(CA»ga.)->CA*E 

AB«IM  X 

^^<^HCAOt .  REAP  1  ^ 

RMidl  1  ••A*«r  r«£«r^ 
(rovn  lO 

«TACP  +  a 

Xw  W«5)-»9TAC3>+S 
VW  .^STACP^V 


(NCAKP) 


Write  oo 
flenowritert 

•  euB  or  STATION 

data  carps' 


_ _ 

Ac  <4^)-*XLAMeA 

SINTM 

C-^CAPC  S— CAPS 

H  <Asi^!ii"T^AL.T 

z  — »  CA-PM 
C0»4  -*■  C0SPH 


Test  ^ 

^NTAPCC) 


RChiSrWl  *F!xpA  6 
I  C1SAVC6) 

MTAPe* 


RcmI  t-U  ,  -C 
6rMn  Dape  6  int* 
TXY*  M<ruTXyf  4-3 
arid  Ji  in+o 
^'S9C3>  fkriA 

(  NTAPC6)— 1 
-*•  NTAPC6 


(Ac  +  ®|tw)  •'*  ««* 

-•  XCMPTH 
eCr¥W)-»THTA 
CMe  •—  C^STH 
X  -♦  CAW 
Oh.  HH0X 
SlV«e-«>6tKTM 
Y  CAFY 
/)-  -P  RH0y 
oj  -p  RM0E 
yO  RH01 
Uk'*><m,>U.Y.XUT 


^L>K<0\ 

Cl.C.  I»  Miellltl^ 

w  U&IaiiU 


A  W«j)  -*•  CAPA 
tiC4lS>-*  SMU4 

doM  fiO-^CMTWPrO 
mint.  Ia  Jnu-iyC^JTWPrrt 
h<\0  •  -P  trCMP 
4a^  ITCMPZ 
Ar  -p  rrcMpa 

•Wi'rx .  -♦  tTtlAP4 
lOOOO  MCrrTDtPZ 
txCloOPP**'  of  Aioa) 

— pTCMP 

A  Cio^ooHk  »f  A  4m) 

I  -PTCMF3 
pC<AAi>tn)-*TCtn>4 


U  -4irXL.SOBX, 
KL%UI^Y,XLSUBZ 

SCJ^) ®euta 

«  (Krs-)-*  AL-PM  A 
^  — r  CPXW 
S'  -P  CFY3>T 

P  (•P4il/|t£'««A) 

-»  ■RRATC 


^-Tact 

(line) 


lao 

^  '  '  ■'"  S 

HEAP 
PriAt  p«6« 
k«Adll*^S 


^  QL0P.  CL®^  ^ 
Wrif*  “■'»  ■IAf«5‘  *■ 

««•«««•  AO.. A. 


l^^TeetN 

Ottiir) 


lPur«:.\A.P«Mfr  — => 

I  GUOP.  eNO  ■*-VH 

I  ^  -  ->  GUOP.tS  ^ 
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BocV^W  \ 

_-Aert  louse  01 


CtV.csm.-kAW^ov  Utt4 


S.O  -^tpfg 

p\us  t  <^*<44  0+ 
^OLrba&e  o>v.tatfC.  ( 
getoiyg  ^e.  (, 

Pot  ttve  corrt.  4Tv<n 

,  buiftr  vuJa  of^refriojTe 

IpAi  l-n.  DLy«*\  Hvru. 

t)i-Tf«+s-  iP«V^. 

torr.  tu±  v-»v  f/q 

|J  * 

c^e  soM-w^  sofAa- 
KvA/sec  oeLTft-h/R 
oeuft-tvw  oeeTft-KYN 
OEltr-uo  ocvTft-Atooe 
OEv-TR-r"  — •>  cKfet> 
tV\rvJL  C.«feO-+n 

HEPiD  N 
ifVivt  fxx^e  j 
\Vv«od^u^_y 

Up&cut')-»  V  — »xvec>^^ 

(.  - •>  UXNfc 

~ 

V  <J«£.K)T)  — 

V  tecOTa) 


PVi-vA:  ccrreetvcM^ 

^  to  tV\e,  cXe-Mwc.'^^ 


Isuka') 

■7K  =(^«TtVL 


CpMCOT) 


s=  COMC-O.') 


CpMtWT  )  -*0MCTa. 
p/o  — •»  OKCWT 

— >G.«TfeV 
UyMa)<jivsTSfc) 
— » C.wr^a. 


fl\.(.3y»<'^ 
v^Test 
(.GnTta)  ^ 

CVW3 


(Sowa^  -^feflTfca 


tfia 


rto  — *  eF«.«6 

I  ttTXHfeS)  -V 
I _  — •  xrtMts 


Vest  ^ 
AtTTMfcS') 


CuftCO 

t  Covnpute.  'v\e>o  sits  t 
-*KXW 

I  ^ 

^  -  tC>l.TMO^vU«6Of+0 
,  — ^  'tnooe 

X  — ►XtblCL 

Cot  x'  —>>  tost 
wi  — ♦vwt 

•lo  —  *1-0 

&  Ekstovv 
«.iw.x'  — >  <StMt 
Cos  So'  —o  c.ose 
w&'  — ^^iKo 

VOy'  — bMOK  wx,'  -»V»W 

e‘*--^es.Q. 

y  — »P,Y^— *KrP. 

Vvy  - *■  rtt  eiHXt 

VvL  — b  Htvi  tMurvt 
xi:  — vHtvalHtrA-a 
Q? 


^  VNISB  S 
Co-vnfOtC  t  ,t. 


(jtOUCvJ.')  vt  ,Tl.*l 
(jynStl  VV.  .TIV 

— =^pa&£_ 

('^"PRtvjrt  ^ 

?r\-»<X.  ■WtVO  1 

oXe-'c^^  J 

Ts^UZZIl 

Utl^TStVx'i  -» 
x\s.-ret+X- V 


■'vest^ 

LtKPT) 


n4-  I 

CWST^-i-oytRPT 
VuroA  e^f  att 

♦  ^oxk.  fer 
/'^\-«Aeict 
fS€C3)e««nr. 


SL'*Vs. 

tVere  <x  X^ 
^CstoLft 


6lr^Uet>-»T>St 

OutyA  vet\A>t 
V'xWTwotienA. 
ovxt'o  'tape  S 


r^OUTINC  TO  KCAO  OBtCKVATION  CARDS 


nAe  \ 


Set  tK\t  Rcnwix 
SoNt  *>1  '•«.  «.3 

Xj 

V».X.Vi  ori^ivioJU 
=xajFvbe^' 


.  r)C.R .  OCR 
f  ^Vo  CO W .  t- 
0+  C-b-s  coorri 
^Cav  -tX^t  v^\.  ^ 
^(6V‘I>  -» 


ec  ^-b  ‘o^io<.ee^4k) 
■»•  ►SRTtcl 

11 -1^  <JBOO©X>J<xy 
-  bt^vbo'-  ro  -»ST<\XO 
->XX1. 

MO-aCK — »  MOf^TH 

- b  OD 

VxC - ?>  HH 

W'nS.  — »  XMM 
bec&.  -  *  '&I 


ra^iA^ 

Pt'cXvO  — >  «?.Rmge 

C  D1  +  L 
^lo  — ^►T^cPhih^ 

0e..Tft,R(5O-q5 
t’t  — * Dl 


''cctTotR.ia 
^Pa-t’iAPf)  — 
R.«>^G£  %/  DL+t 
A  bCcs.tb  tiwe 
-»&L 


<iex 

>-* 

WrtcJk 

ti  3 

octf.  OC*V3> 

e.A)iG>E  2<  OL-vl 
•ViA.  a  sect-t)  f IV 

^  R  cm 


CX«.C>CRi<l 
r  rvocebt  d^b  cr  ft,V 
— ArvR 

^SrrS.— >0ELTR  ^ 

^<0 

<b<.<.oce3 

'Vnstebb  1.  Ooff>\e.r  ^ 

.  T/c  — ^  RuPHP*, 
JbEt-TH.RftKlC^ 

Rt-RPA- 


r„  cce.ex,Ri.v 
ri.=  ).\ireA 
«.ftKNc,-.E  S  Pt+l 
Hdd  A.  tr  twAC 

t#>v  — *  PI. 

ft  ovc  exi-fc 


n»»t.pb 

\0p'^  oV 
--i.CFRgCXV 


,  OCK.OCRlt 

R-RNSE  t  OL+1 

ftAck  a  veci.t;tmey 

-V  Pit 
ft  caACxVt^ 


ote.ociRvs’ 

Pa 

R.R>4G6  «>  OJ.-*! 
VVid  a  secaU-ti-VAe 
^PtU.«*0  Dt 
V  ft  O-A  <ix<t^ 


vN» 


[ _ 

|Sci  %vjO\'\dh  to  1 

m 

■^c.  ocRn 

ftftNGE  <2,  Pt  +  \. 

ftAA  a  \ec.s. 

TVvv\g 


u.  'more  Cords 
Vrowv  e  O  \-v\So 
XS'JF  YVu.  XeuF^  \tft 
X60F  -H»X. 


R.CRRDX 


, 

RfioA  \a  >«or6  Carris 
•Wpvn.  tn^ie  0  VAiCto 
X60F  «rru  X0JF+Vtft 

r  r  i 

X.B0F  — »-X. 

l.X,l— ^K1 
■L  -♦  NePlRD 

, 

tR-estbr-e  X., 

- 1 _ 

^rp-r*\  IIC3 

M 

A  — »  MCftftO 

.44. 


ROUTINC  TX>  PMCf^APUl  OB»CKVATION  RECORD 


Compute : 

p— U-P 

CSQ. 

a— pA  JeT KTA 

— •  XN  — •  'KTP 

U  — »  Wx  ^ru  WX-t-2 

co*£.  — ♦  ctfsl 

tint  — •  «|NZ 
M^-*’KMZ  M,-»XNX 
Ny-rXNy 

Mh-^XMX 
Clx"  —  AXNI 

-►  ayni 

XN0OC 
U  <r*a)— •  U 


vcn»a) 


U(>Wiar)-»  e^2 


f(B^-(0i)r 

.  ^ .ooooolv 


^C+w' 

InCIPt 


•  irtCC'Kij)-^  SiNC# 
e»>  Cc.-KJ)-»  C^E# 
Cftin  C-^  ESiHC 
ECOtC-PCC^SC 
r-—  w  Vr  -•  AX 
ft-*.*  — p  XTCSQ. 

I  +  Zirs?-^  XLSR 
»'A  u.  — »  CIMU 
cepu  — ►  C0SU 

J^-^UK  Ikrut.  UX*2 
^-pVX  VXV2 

-p  ■♦X  Hr«*  X'fZ 
r  -p  KO0T 
«*'>  -P  RV^T 

■fc  XS#r-«>£ 


Cc^c)~»g^t 


>0  iiV 


Writ*,  en 
Pi*i)*iiijrif«r : 

*«0  TINES  tnhousn 
LOOP  MiTMOUr 
CLOSINS  ONE* 


46- 


RouriNe  TO  comrutc  one  row  or  coefficient* 

FOR  LEAST  SQUARES  MATRIX 


CORRP 


t  dotted  I 

i-«t«  u,v,y  ,a.>ul  ^ 

— ftPOTU  ,  hOoTv, 
ftOOTw .  ftpoei. 


A  — *  RHOOT«.l 
tVvrOL  ^HObT-^a 
RHOOT 
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2.4  INPUT  AND  OUTPUT  FORMATS 

The  Orbit  Correction  Program  uses  two  input  information  tapes, 
Control  and  Observation  Cards  on  logical  tape  unit  0  and  Station  Informa¬ 
tion  on  logical  tape  unit  10.  The  output  is  available  as  printed  output 
of  various  types  on  logical  tape  unit  5,  the  residuals  can  be  punched  from 
the  tape  on  logical  unit  5  with  data  select  2.  The  tape  on  logical  tape 
unit  6  contains  the  positions  and  velocities  at  all  observation  times  for 
use  in  the  Station  Locator  Program. 

2.4.1  INPUT 

Before  transferring  control  to  this  program,  the  data  cards 
described  below  must  be  put  on  logical  tapes  0  and  10.  This  is  usually 
done  by  appropriate  control  cards  preceding  these  data  cards  in  the  deck 
submitted  at  execution  time. 


Control  and  Observation  Cards  on  Logical  Tape  0: 


Four  cards  for  beginning  the  ephemeris  computation  are  needed  plus  1 
control  card  to  specify  the  use  to  be  made  of  the  program.  These  five 
cards  are  all  that  are  necessary  for  Simulation  but  if  Differential  Cor¬ 
rection  is  desired,  these  five  cards  are  followed  by  any  number  of  obser¬ 
vation  cards  which  are  followed  by  an  END  card  (a  card  with  "ENDA<2kAAA" 
in  columns  1-8).  The  formats  of  these  cards  are  described  below. 


Control  Cards  : 

Card  No .  Card  Cols.  Contents 

1  1  Print  option: 

0  means  print  off-line  none  of  the  follow¬ 
ing. 

1  means  print  off-line  t-t  (time  since 
epoch)  in  minutes,  <f)  (latitude)  in  degrees, 
Ag(East  longitude)  in  degrees,  and  H 
(height  above  earth)  in  meters  for  each 
point  of  the  ephemeris. 

2  means  print  off-line  t-t  in  minutes,  x, 
y,  z,  X,  y,  z  for  each  point  of  the  ephem¬ 
eris.  X,  y,  z  are  in  earth  radii  and  x, 

* 

means  blank 
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Card  No. 
1 


Card  Cols. 


Contents 


37-48 

49-60 

61-80 


-1 


y,  z  are  in  earth  radii  per  min. 

means  print  off-line  the  output  for 
both  1  and  2. 


2-6 

A  CASE 

7-12 

AOXXXA --  where  OXXX  is  the  case 
number . 

13-72 

First  line  of  heading  for  each  page 
of  output. 

73-80 

Not  used. 

2 

1-72 

Second  line  of  heading  for  each  page 
of  output. 

73-80 

Not  used. 

3 

1-12 

^400000000+00 

13-24 

At  in  minutes--time  interval  used  in 
Runge  Kutta  integration  of  ephemeris 

25-36 

l  t,-t  in  minures--time  interval  betwef 
\  f  o 

floating 
point 
(see  note2) 


epoch  time  and  final  time  point  in 
ephemeris  integration. 

d  in  meters--the  caliber  or  reference 
diameter  of  the  vehicle. 

m  in  kilograms--the  weight  of  the 
vehicle . 

,Not  used. 


1-12 

13-24 

25-36 

37-48 

49-60 

61-72 

73-80 


in  radians --mean  longitude  at  epoch. 

=  e  cos  CJ 
xN  o  o 

o 


floating  ]a  =  e^  sin  CO 
point  /  ^  o  ' 


Components  in  orbit 
plane  at  epoch  of  a 


(see  note 9  \  h 


Components  of  angular  momentum 
vector  at  epoch,  h^ 


VNot  used . 
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Card  No .  Card  Cols . 


Contents 


5  1-6  XXXXXX--these  characters  determine  which 

parameters  are  to  be  corrected.  Each 
character  may  be  a  "0"  or  a  "1"  where  a 
zero  indicates  that  the  corresponding 
parameter  is  not  to  be  corrected  and  a 
one  indicates  that  correction  is  desired. 
All  zeros  indicate  that  simulation  of  ac- 


quisition  coordinates  is  desired.  The 

order 

of  the  parameters  corresponding  to 

the  6 

1  characters  is:  a  ,  a  „  ,  a  „  ,  U  , 

0  xN  yN  0 

o  o 

n  o» 

i  . 

0 

7-10 

Epoch 

1,  day  of  year--fixed  point  integer. 

11-18 

Time 

from  beginning  of  epoch  day  to  epoch 

in  minutes--fixed  point. 

19 

Not  used. 

20 

Number  of  times  to  repeat  correction. 

21 

0 

22-24 

0"^  (see  note  l)--fixed  point. 

25-28 

Restart,  day  of  year--fixed  point  integer. 

29-37 

Time 

from  beginning  of  restart  day  to  re- 

start  time  in  minutes--fixed  point. 

38-42 

Not  used. 

43-48 

ABSMX  =  absolute  maximum  for  angle  and 

range  residuals--f ixed  point,  (see  note  1) 

49-51 

' 

52-54 

1 

55-57 

58-60 

=^5 

>  fixed  point  (see  note  1) 

61-63 

64-66 

67-69 
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Card  No .  Card  Cols. 


Contents 


5  70-72 


73-78 


■=  Maximum  number  of  times  to  pass 


through  the  observations  before  actually 
correcting  the  parameters.  is  set  to 

3  internally  if  input  as  0.  (see  note  1) 
fixed  point. 


ABMX2  =  absolute  maximum  for  range  rate 
residuals--f ixed  point.  (see  note  1) 


79-80 


Not  used. 


Note  1 


During  differential  correction  there  is  an  automatic  feature  for 
rejection  of  bad  observations  which  works  as  follows:  usually  at  least 
two  passes  are  made  through  all  the  observations  before  actually  correcting 
the  orbit  parameters.  During  the  first  pass,  before  correcting  the  para¬ 
meters  the  i*-^  time,  the  residuals  are  compared  with  ABSMX  for  range  and 
for  angle  observations  and  ABMX2  for  range  rate  observations.  All  the  re¬ 
siduals  whose  absolute  values  are  less  than  these  limits  are  used  in  the 

N 

formation  of  two  root -mean-square  values:  rms^ 


where  the  r , 


are  the  N  range  and  angle  residuals  that  pass  the  test,  and  rmS2 


where  the  s  are  the  M  range  rate  residuals  that  pass  the  test.  During 
the  second  and  subsequent  passes  through  the  observations,  before  correct¬ 
ing  the  parameters  the  i^^  time,  the  range  and  angle  residuals  are  compared 
with  O'  times  the  rms  of  the  previous  pass,  while  the  range  rate  residuals 
are  compared  with  O’ ^  times  the  rms2  of  the  previous  pass.  As  before,  only 
the  residuals  that  pass  the  test  are  used  in  computing  the  next  two  rms 
values.  Only  the  observations  whose  residuals  pass  the  test  are  used  in 
the  computation  of  the  corrections  to  the  orbital  parameters.  The  para¬ 
meters  are  not  actually  corrected  until  either  the  same  number  of  residuals 
are  rejected  on  successive  passes,  or  the  maximum  allowed  number  of  passes 
(  O’  )  made . 
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Note  2 


By  a  "floating  point"  number  is  meant  a  number  in  the  following 

format : 


txXXXXXXXlYY  (this  is  for  a  12  column  field) 

where  ^XXXXXXXX  is  the  mantissa  and  IVY  is  the  power  of  10  that  it  is  mul¬ 
tiplied  by.  The  decimal  point  is  assumed  just  before  the  first  X,  which 
must  be  non  zero  (unless,  of  course,  the  number  equals  zero).  For  instance, 
-25  would  be  written  as  -25000000+02  and  +.0003  would  be  written  as 
+30000000-03  or  30000000-03. 


Observation  Cards: 

These  cards  must  be  ordered  with  respect  to  time.  There  are  four 
different  types  of  observation  cards:  Radar,  Optical,  Doppler,  and  Emitted 
Frequency  cards.  The  information  on  the  first  three  of  these  types  of  cards 
in  columns  1-29,  is  the  same  and  is  as  follows: 


Card  Cols. 


Contents 


1-2 

3-4 

5 

6 

7-10 

11-14 

15 

lG-17 

18-19 

20-21 

22-23 

24-25 

26-29 


Year  of  launch  (1950=01) 

Satell:ite  Greek  letter  (  0<  =01,  etc.)  Satellite  identifica- 
Satellite  component  number  tion 

Observation  type  (1  =  Radar,  2  =  Optical,  3  =  Doppler) 
Observation  or  pass  number  (not  used  by  program) 

Station  identification 
Year  (1958=1)  \ 

Month 

Day 

Hour  \  Observation  time 

Minute 

Second 

.0001  second  y 
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The  observation  time  must  be  referred  to  the  besinnina  of  the 
same  year  as  the  epoch.  Days  in  January  of  the  next  year  should  be  re¬ 
ferred  to  as  days  of  December;  1962  January  15  =  1961  December  46.  0 

S’^o 

(Greenwich  Sidereal  time  at  epoch)  is  computed  within  the  program.  To 
do  this,  Greenwich  Sidereal  time  at  the  beginning  of  January  0  of  epoch 
year  is  a  constant  in  the  program.  At  the  time  of  this  writing,  tliis  con¬ 
stant  is  for  January  0  of  1961  (99?420937).  Therefore,  to  use  an  epoch 
that  is  in  any  other  year  than  1961,  requires  a  change  of  this  constant. 


Radar  Observation  Card: 
Card  Cols. 


1-29 

30-34 

35-42 

43-47 

48-52 

53-57 

58-62 

63-67 

68-72 

73-80 


Contents 


Identification  and  time  of  observation. 

Year  and  day  of  reduction  (not  used  by  program) 

Slant  ran.ge  in  meters--f ixed  point  integer 
Difference  in  slant  range  from  t  to  t+2  seconds 
Difference  in  slant  range  from  t+2  to  t+4  seconds 

Difference  in  slant  range  from  t+4  to  t+6  seconds  I 

Difference  in  slant  range  from  t+6  to  t+8  seconds  /point 

Difference  in  slant  range  from  t+8  to  t+10  seconds\ 

Difference  in  slant  range  from  t+10  to  t+12  seconds) 
Slant  range  at  t+12  seconds  (not  used  by  program).^ 


Optical  Observation  Card; 
Card  Cols. 


1-29 

30-32 


33-34 

35-36 


Contents 


Identification  and  time  of  observation 

Right  ascension  (hours)  or  azimuth  (degrees) - -negative 
sign  in  column  30  means  that  this  is  a  right  ascension 
and  declination  observation,  anything  else  in  column  30 
means  an  azimuth  and  elevation  angle  observation. 

Right  ascension  or  azimuth  (minutes  of  time  and  arc, 
respectively) . 

Right  ascension  or  azimuth  (seconds  of  time  and  arc, 
respectively) . 
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Card  Cols. 
37-38 

39-42 

43-44 

45-46 

47 

48-58 

59-80 


Contents 

Right  ascension  or  azimuth  (.01  seconds  of  time  and  arc, 
respectively) . 

Declination  or  elevation  (degree)  (right  adjusted). 
Declination  or  elevation  (minutes). 

Declination  or  elevation  (seconds). 

Declination  or  elevation  (0.1  seconds). 

Miscellaneous  information  not  used  by  program. 

Not  used. 


Doppler  (Range  Rate)  Observation  Card; 


Card  Cols. 


Contents 


1-29 


Identification  and  time  of  observation. 


30-37 

38-48 


49-53 


Miscellaneous  information  not  used  by  program. 

Received  frequency  in  hundredths  of  cycles  per  second, 

scaled  to  108  megacycles  per  second,  i.e.,  108,000,000  must 
be  added  to  obtain  actual  received  frequency  in  cycles  per 
second  (i.e.,  this  field  contains  )/*  -  108  me.). 

Miscellaneous  information  not  used  by  program. 


54-80  Not  used. 

• 

To  compute  range  rage.  A,  from  the  above  information,  the  follow¬ 
ing  formula  is  used:  p  =  -c  -  where  c  is  the  speed  of  light  and  is 

the  frequency  emitted  By  the  satellite.  J/is  obtained  from  the  most  recently 
encountered  Emitted  Frequency  Card  (format  described  next). 


Emitted  Frequency  Card 


Card  Cols. 


1-5 


Not  used. 


Contents 


6 


Blank--this  differentiates  this  card  from  the  other  obser¬ 
vation  cards. 
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Card  Cols . 


Contents 


7-29 

30-40 


41-80 


Not  used. 

Frequency  emitted  by  the  satellite,  ^  ,  in  hundredths 
of  cycles  per  second  scaled  to  108  megacycles  per  second 
(l.e.,  this  field  contains  V  -  108  me.) 

Not  used. 


b .  Station  Information  Stored  on  Logical  Tape  10: 


Any  number  of  station  cards  in  any  order  followed  by  an  END  card 
(a  card  with  "ENdAAAAA"  columns  1-8),  can  be  read  onto  tape: 


Card  Cols. 


1-4 

5-12 

13-21 


22-27 

28-80 


Contents 


Station  identification. 

0,  latitude  in  one  hundred  thousandths  of  a  degree,  l.e., 
there  is  a  decimal  point  assumed  between  columns  7  and  8. 

Xg,  east  longitude  in  one  hundred  thousandths  of  a  degree, 
l.e.,  there  is  a  decimal  point  assumed  between  columns  16 
and  17. 

H,  height  above  sea  level  in  meters. 

Not  used. 


For  simulation,  acquisition  coordinates  will  be  computed  at  all 
points  of  the  ephemeris,  for  all  stations  which  can  see  the  vehicle. 

During  differential  correction,  if  cols.  11-14  of  an  observation 
card  do  not  match  exactly  cols.  1-4  of  some  station  card,  the  observation 
will  be  skipped  and  "STATION  NUMBER  OOOOXXXX  NOT  FOUND"  will  be  written  on 
the  flexowriter. 


2.4.2  OUTPUT 


The  output  of  the  Orbit  Correction  Program  is  used  to  monitor  the 
correction  and  simulation  procedures,  to  provide  restart  information  when 
the  correction  procedure  is  interrupted,  and  to  provide  input  to  the  Station 
Locator  Program.  The  following  describe  where  each  type  of  information  is 
to  be  found. 
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a .  Printed  Information: 


The  following  printed  output  is  obtained  by  printing  tapes  with 
data  select  1.  This  output  depends  considerably  on  whether  simulation  or 
differential  correction  is  being  done.  However,  there  is  some  output  that 
is  the  same  in  either  case.  All  pages  are  headed  by  the  information  on 
cards  1  and  2  (see  INPUT).  If  the  print  option  (col.  1,  card  1)  is  equal 
to  1  or  3 ,  t-t  (time  since  epoch)  in  minutes,  (latitude)  in  degrees, 

Ag  (East  longitude)  in  degrees,  and  H  (height  above  earth)  in  meters  are 
printed  under  appropriate  headings  for  all  points  of  the  ephemeris.  If 
the  print  option  is  equal  to  2  or  3,  t-t 
(x,  y,  2  in  earth  radii  and  x,  y 


in  minutes,  x,  y,  z,  x,  y,  and  z 
2  in  earth  radii  per  '' 


^min)  are  printed 


under  appropriate  headings  for  all  points  of  the  ephemeris. 


After  Simulation: 

In  addition  to  the  above,  t-t^  in  minutes, (range)  in  earth 
radii,  p  (range  rate)  in  earth  radii  per  k  min,'cx  (right  ascension) 
in  degrees,  q  (declination)  in  degrees,  A  tazimoth)  in  degrees,  and  h 
(elevation  angle)  in  degrees  are  printed  under  appropriate  headings  for 
every  station  for  points  of  the  ephemeris.  Under  each  station  only  those 
points  of  the  ephemeris  are  output  where  h  >  0  (hethe  satellite  is  above 
the  horizon). 


After  Differential  Correction: 

For  all  passes  through  the  observations  all  residuals  are  printed. 
Printed  with  each  residual  are  the  satellite  identification  (cols.  3-5  of 
the  observation  card),  the  station  number  (cols.  12-14  of  the  observation 
card),  and  the  time  since  epoch  in  minutes.  All  quantities  are  listed 
under  an  appropriate  heading.  The  range  and  angle  residuals  are  in  kilo¬ 
meters,  while  the  range-rate  residuals  are  in  kilometers  per  second.  The 
rejected  residuals  are  flagged  at  the  right  of  the  line  by  four  asterisks 

(***yf)  _ 

On  the  next  page  following  the  residual  information  for  each  pass 
are  printed  the  case  number  (cols.  8-11  of  ephemeris  card  1),  the  root-mean- 
square  of  the  good  range  and  angle  residuals  in  kilometers  (under  the  head¬ 
ing  SUM-KM) ,  the  root -mean-square  of  the  good  range-rate  residuals  in  l.ilo- 
meters  per  second  (under  the  heading  SUM2-KM/SEC) ,  and  the  corrections  to 

the  orbital  parameters,  Aa  /  ,  Aa.,  ,Aa„,  AU  ,  A  >  and  .A  i  , 

o/a  xN  yN  ’  o  o  o 

based  on  the  good  observations  from  this  pass. 
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After  the  final  pass  through  the  observations  when  these  correc¬ 
tions  are  actually  added  to  the  orbital  parameters,  the  following  informa¬ 
tion  is  printed  just  after  these  corrections.  The  case  number  (which  is 
increased  by  one  in  the  left-most  position  at  each  correction)  and  the  new 
orbital  parameters  at  epoch,  Lq,  a^N  ,  ''yo’ 

y,  z.  y>  z  at  epoch  are  printed. 

After  all  the  corrections  have  been  done,  if  there  is  a  restart 
time  specified  in  the  control  card,  then  the  orbital  parameters  at  this 
restart  time  are  printed. 

b.  Punched  Information: 


The  punched  output  is  obtained  by  punching  tape  11  with  data  select 
2.  There  is  no  punched  output  for  simulation.  In  the  differential  correc¬ 
tion,  all  the  good  residuals  from  the  last  correction  are  punched  and  may  be 
used  as  input  to  the  Station  Locator  Program.  The  format  is  described  in  the 
Station  Locator  Program  input  specifications  (Section  3.4). 

On  Binary  Tape: 

At  the  end  of  differential  correction,  the  tape  on  logical  tape 
unit  6  contains  t-t^  in  minutes,  x,  y,  z,  x,  y,  and  z  at  all  the  observa¬ 
tion  times  (x,  y,  z  are  in  earth  radii  and  x,  z  are  in  earth  radii  per 

k*  min).  This  tape  is  to  be  used  as  input,  in  conjunction  with  the  above 
residual  cards,  to  the  Station  Locator  Program, 


2.5  OPERATING  PROCEDURE: 

The  deck  submitted  at  execution  time  has  the  following  make  up: 

Starting  in 

Col.  17  Col.  25  Purpose 

JOB  XXXXX...  Signals  the  start  of  a  new  job.  Cols.  25 

on,  are  typed  on  the  flexowriter  at  the 
beginning  of  the  program  and  may  contain 
any  alphanumeric  information. 
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Col.  17 

Col.  25 

Purpose 

REWIND 

0,3,6,9,10 

Rewinds  the  indicated  tapes. 

RPL 

1 , DATA, GO 

Reads  routine  named  DATA  from  tape  1 
memory  and  transfers  control  to  it. 

into 

TAPE 

0 

Control  instruction  for  DATA, telling 
which  tape  to  put  the  following  data 

it 

on . 

Data  to  be  put  on  tape  0  (cf.  Sec.  2.4.1) 


ENDDATA 

Control  instruction  for  DATA  telling  it 
that  there  is  no  more  data  to  be  put  on 
tape  0. 

REWIND 

0 

Rewinds  tape  0. 

JMP 

* 

Transfers  control  to  the  program  in  memory 
(DATA) .  (See  note) 

TAPE 

10 

Control  instruction  for  DATA  telling  it 
which  tape  to  put  the  following  data  on. 

.... 

Data  to  be  put  on  tape  10  (cf.  Sec.  2.4.1) 

ENDDATA 

J 

Control  instruction  for  DATA  telling  it 
that  there  is  no  more  data  to  put  on  tape  10 

REWIND 

10 

Rewinds  tape  10. 

RPL 

4, XXX. 

.  .XXX, GO 

XXX... XXX  is  the  identity  of  the  Orbit  Cor- 

rection  Program  (up  to  16  characters)  .  It 
is  read  into  memory  from  tape  4  and  control 
is  transferred  to  it. 
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Note  : 


The  card  pertaining  to  tape  10  may  be  omitted  during  differential 
correction  If  a  previously  prepared  station  tape  has  been  mounted  on  logical 
tape  10. 


Tapes  Used  : 
Logical  Tape  No. 
0 

3 

4 

5 

6 


9 

10 

11 


Purpose 


Data  Input. 

Ephemeris  tape--used  to  save  t-t^,  L,  a, 
h,  at  all  ephemeris  points. 

RPL  tape  containing  Orbit  Correction  Pro¬ 
gram. 

Output  tape. 

Used  during  computation  of  ephemeris,  during 
simulation  or  during  differential  correction, 
if  the  print  option  is  2  or  3,  to  save  t-tg, 
r,  r  at  observation  times  to  be  used  by 
Station  Locator  Program. 

Observation  tape  --observations  are  processed 
and  saved  on  this  tape. 

Data  input.  (Station  tape) 

Last  correction  residuals 


To  Use  the  Orbit  Correction  Program  with  the  Station  Locator  Program: 


(1)  Run  the  Orbit  Correction  Program  with  nominal  station 
coordinates . 

(2)  Save  tape  6--binary  tape  containing  t-tg,  r,  and  £ 
at  observation  times. 

(3)  Print  tape  5  off-line,  using  DSl  (data  select  1)  for 
printed  output  from  Orbit  Correction  Program. 
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(4)  Punch  tape  5  off-line  using  DS2  (data  select  2)  for  the 
general  information  card  (Ogro^  input  to  the 
Station  Locator  Program. 

(5)  Punch  tape  11  off-line,  using  DS2  (data  select  2)  for 
punched  output  from  Orbit  Correction  Program.  The 
residuals  are  punched  for  the  last  pass  through  the 
observations. 

(6)  If  desired,  punch  (by  hand)  weights  in  the  residual 
cards.  If  none  are  punched,  I's  will  be  assumed  by 
the  Station  Locator  Program. 

(7)  Punch  additional  input  cards  for  Station  Locator  Pro¬ 
gram  (station  information  card(s)  ),  and  set  up  input 
deck. 

(8)  Run  the  Station  Locator  Program  with  the  binary  tape 
containing  t-tg,  £  and  r  at  observation  times  prepared 
by  the  Orbit  Correction  Program  on  logical  tape  3. 

(9)  Print  tape  5  off-line,  using  DSO  (data  select  0)  for 
printed  output  from  Station  Locator  Program. 

(10)  If  desired,  repeat  Orbit  Correction  with  new  station 
coordinates,  to  obtain  new  elements  and  new  residuals. 

(11)  Then,  if  desired,  repeat  Station  Locator,  to  refine 
station  corrections. 
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SECTION  3 

STATION  LOCATOR  PROGRAM 


The  purpose  of  the  Station  Locator  Program  is  to  compute  from 
observed  deviations  by  a  geodetic  satellite  from  a  good  model  of  its  orbit, 
corrections  to  the  geocentric  coordinates  of  the  stations  observing  the 
satellite  and  to  the  geocentric  coordinates  of  the  origin  of  any  datum  to 
which  a  group  of  such  stations  are  tied.  The  good  orbital  model  results 
from  the  Orbit  Correction  Program  and  the  observed  deviations  are  the  resid¬ 
uals  produced  by  that  program. 

Since  the  station  coordinates  have  been  determined  already,  with 
a  certain  accuracy  from  previous  information  (e.g.,  surveys),  it  is  not 
desirable  to  base  the  determination  of  the  corrections  entirely  on  the  ob¬ 
servations  of  one  geodetic  satellite.  Therefore,  limits  are  placed  upon 
the  magnitude  of  the  station  corrections.  Furthermore,  there  is  provision 
for  weighting  the  stations  on  the  same  datum  to  account  for  different  strengths 
of  their  ties  to  the  datum  origin.  The  individual  observational  quantities 
may  be  weighted  at  will,  to  account  for  observing  conditions,  etc. 


3 . 1  THEORY 

Assuming  the  orbit  of  an  earth  satellite  is  known,  any  observation 
of  the  satellite  can  be  expressed  in  terms  of  the  time  and  three  coordinates 
describing  the  location  of  the  observer.  Equations  (1),  (2),  and  (3)  of 
Section  2.1,  are  again  applicable  where  the  X.  are  now  the  three  station  co¬ 
ordinates  and  (C.  .)  is  an  m  X  3  matrix.  The  '^elements  of  this  matrix  are 
1  J 

now  the  partial  derivatives  of  the  observed  quantities  with  respect  to  the 
station  coordinates.  The  solution  again  takes  the  form  of  equation  (4),  and 
thus  the  iterative  process  would  result  in  improvements  to  the  station  coor¬ 
dinates. 
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When  weighting  is  desired  in  the  determination  to  account  for 
different  types  of  tracking  data  or  different  tracking  instruments,  the 
least  squares  matrix  can  be  constructed  according  to 


"ip> 


(6) 


where  (P.  ) 
iP 


is  the  diagonal  weight  matrix,  defined  as 


P. 

J-P 


=  . 
oi 


i  =  p 


P.  =  0  i  jt  p 

i-P 


The  cO  .is  the  weighting  factor  of  the  observation  0. .  The  solution  now 
tah.es  ?fte  form:  ^ 


(^X.)  =  N'^  (0^.)“^  (P^p)  (Ao^) 


(7) 


The  A  X  are  compared  against  limits,  chosen  with  regard  to  the 
uncertainty  in  that  coordinate  prior  to  the  geodetic  satellite  observations. 
If  one  of  the  corrections,  say  AX^,  exceeds  its  limit,  the  correction  is 
set  equal  to  the  limit  with  the  sign  of  the  discarded  excessive  correction. 
The  result  is  to  transform  the  equations  of  condition  from 


,>o  0 

Aoj  ■  Ax. 


to 


Ax 


A 


(8) 


(9) 


This  set  of  equations  can  be  solved  for  A  X^  and  AX  as  previously  shown, 
where  the  (C^,)  matrix  is  now  an  m  x  2  matrix,  and  the  least  squares  matrix 
N  becomes  a  2x2  matrix. 
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In  the  . teratlve  process,  the  corrections  to  the  two-station 
coordinates,  as  obtained  from  equation  (9),  would  be  compared  to  their 
limits.  This  test  determines  whether  the  solution  is  acceptable  or  whether 
a  "least  squares"  solution  in  terms  of  only  one  of  the  station  coordinates 
is  more  desirable.  Such  a  one  coordinate  solution  amounts  to  a  weighted 
average . 


The  following  example  shows  a  fictitious  case.  The  units  are 
given  in  the  last  column. 


Solutions 


Limits 

1 

2 

3 

A  X^ 

= 

+  0.003 

+  0.00177 

H  0.00279 

-  -  degrees 

■  ^2 

=  cos  1  AA 

+  0.003 

+  0.00045 

H  0.00192 

-  -  degrees 

A.X3 

=  A.H 

t  100 

-  327 

-  100 

-  -  meters 

A  Xj  from  the  first  solution  exceeds  its  limit.  It  is  set  equal 
to  that  limit  with  the  appropriate  sign.  Since  neither  nor  AX  ,  from 

tl'.e  second  solution  exceed  their  limits,  the  third  solution^is  not  necessary 
and  the  second  solution  is  adopted. 


3.2  LOGICAL  PROGRAM 

The  cartesian  station  coordinates  are  calculated  for  each  station 
under  the  condition  that  the  Greenwich  meridian  is  aligned  with  the  great 

circle  passing  through  the  poles  and  the  vernal  equinox  (i.e.,  ©  =0) 

gr 

X  (L)  =  -  X  cos  A 

s  c 

Y  (L)  =  -  X  sin  A 

s  c 

Z  (L)  =  -  y 

s  c 
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where  y^.  =  (S  +  H)  sin  ^ 
and  *  (C  +  H)  cos  ^ 

The  quantities  C  and  S  are  computed  as  follows: 

-ic 

C  -  ag  (1  -  e^  sin  ^  S  =  (1  -6^)0 

2 

The  constants  ag  (semi-major  axis  of  the  ellipsoid)  and  e  (e  is  the 
eccentricity  of  the  ellipsoid)  are  now  constants  contained  within  the 
program: 


ag  =  1.000  000  000  0 

e^  =  0.006  693  421  6 


A  logical  change  is  necessary  at  this  point  to  allow  for  a  look¬ 
up  of  a^  and  e  according  to  the  datum  number.  It  is  advisable  that  these 
"constants"  correspond  to  the  ellipsoid  used  to  tie  the  datum  together.  For 
each  station  the  constants,  CKl,  CK2,  CK3,  and  CK4  are  calculated  as  follows: 


CKl 

CK2 

CK3 

CK4 


(C  +  H)  +  e^C^  cos^^a"^ 

(S  +  H)  +  e^SC^  sin^()ag^ 

Xg  tan^  -  e^C^  cos^^  sin^ag 

2  2  2  A  ~  2 

yg  cot^  +  SC  e^  sin'^()  cos^^ 


Calculation  of  0,  the  sidereal  time,  from  the  observation  time,  t,  gives 
the  components  of  the  station  location  vector  at  this  time:  XI,  Yl,  and  Zl. 


0  =  0gr  +  0(t  -  tg)  +  > 

Where  0gj.  is  the  angle  between  the  vernal  equinox  the  the  Greenwich 
meridian  at  epoch  time  t^,  ^  is  the  earth-fixed  longitude  (i.e.,  its  reference 
is  the  Greenwich  meridian)  and  d  is  the  rotation  rate  of  the  earth. 


The  station  components  at  this  longitude  are: 


XI  = 

-  ’'c 

cos 

Yl  = 

-  ’'c 

sin 

Zl  = 

-  yc 
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a.  Calculation  of  the  Components  of  the  Topocentrlc  Range  Vector 
The  range  vector,  ^  ,  has  components 

/O  =  X  +  XI 

'  X 

/O  =  y  +  Y1 

.  y 

0  =  2  +  Z1 

•  2 

from  which  the  direction  cosines  are: 

K  ■  P^ip  _ ^ 

Ly  =  where  P  -  ' 

■•e  ■  PJP 

\ 

The  geometric  quantities  x,  y,  2,  the  geocentric  components  of 
the  radius  vector  and  the  satellite,  at  time  t,  were  read  from  a  binary 
tape,  prepared  by  the  Orbit  Correction  Program. 

b .  Calculation  of  the  Weighted  Differential  Coefficients  for  the 
Different  Types  of  Residuals 

The  residual  in  range,  AP,  is  related  to  errors  in  ^  ^  , 

and  H  as  follows:  ' 

Ap=R^A<^  +  R^  cos  A'X  +  R^  Ah  (10) 

Where 

R1  =  (L  cos  ©  +  L  sin  ©• )  CK3  -  (L  )  CK4 
X  y  2 

R2  =  (C  +  H)  (L^  sin  ©  -  cos  ©  ) 

R3  =  -  cos  d)  (L  cos  ©+L  sin©  +L  tan  <b  ) 

T  X  y  2  ' 

At  this  point  the  weighting  factor  must  be  taken  into  account,  and  equation 
(10)  becomes 
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The  coefficients  in  this  weighted  equation  are  the  quantities  CE(1)  ,  CE(2)  , 
CE(3)  and  RES,  defined  as: 


CE 


(1)  =  R^ 

CE(2)  =  R2 

CE(3)  =  Yv^  R3 

RES 

These  quantities  are  then  utilized  to  construct  the  least  squares  matrix. 
This  is  done  in  the  portion  of  the  program  referred  to  as  subroutine  Sub  1. 
Utilizing  the  same  theory  as  stated  above,  the  following  relationships  are 
obtained  for  the  residuals  ^  ,  PAA  cos  h,  /®Ah  and  A/®. 

The  quantities  CE(1),  CE(2) ,  CE(3)  and  RES  are  calculated  in  precisely  the 
same  manner  as  they  were  calculated  for  the  residual  of  AP,  except  that 
the  differential  coefficients  depend  upon  the  type  of  residual  being  uti¬ 
lized,  and  will  thus  vary  from  equation  (10). 

For  the  residual  ^A  \  cos  ^  : 

/O  AoC  cos  S  =  (aLI  Af|^+  AL2  cosC^)  +  AL3  AH^ 

CE(1)  =  AL) 

CE(2)  =  Yw^'  AL2 

CE(3)  =  AL3 

- 


or 


where 


RES  =  /  ^A«Kcos 
ALI  =  CK3 


CA  cos  ©  +  A  sin  ©  ) 

X  y  J 


sin  © 


A  cos  ©  ) 

y  J 


AL2  =  (C  +  H)  ^A^  sii 

AL3  =  -  cos  <p  cos  0  +  A^  sin©^ 
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-L 


and : 


A  = 


1  -  L 


T~> 


A  = 

y  1 


1  -  L 


T-i 


where 


and: 


isidual  ^ 

<1 

e 

II 

< 

CE(1)  = 

yr 

DLl 

CE(2)  = 

DL2 

CE(3)  = 

DL3 

RES 

DLl  =  -  CK4  (D  )  +  CK3  (D  cos  0  +  D  sin  0 ) 
z  X  y 

DL2  =  (C  +  H)  (D^  sin  0  -  cos  © ) 

DL3  =  -  cos  (D^  cos  ©  +  D  sin  0  +  tan  <!p  ) 


D  =  -  L  A 
X  z  y 


L  A 
z  X 


1  -  L 


For  the  residual  r^A  cos  h: 

cos  h  =  (azi  a4>  +  AZ2  cos  cj>  A>  +  aZ3  A 

^sin  <(?  (CK3)  +  cos  (CK4)^ 


where:  AZI  =  A 


xh 


AZ2  =  -  (C  +  H)  A 


yh 


AZ3  =  0 
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S^  =  sin  (^5  cos  © 
Sy  =  sin  Ip  sin  © 

=  -  cos 

E  =  -  sin  & 

X 

E  =  cos  © 

y 

=  cos  (^5  cos  © 

Z^  =  cos  <p  sin  © 

Z^  =  sin  (j) 


xh 


JlL 


1  -  L 


*71 

zh 


where 


L 


xh 


L  S 

X  X 


+  L  S 

y  y 


+  L  s 

z  z 


L  E 

X  X 


+  L  E 

y  y 


L 


zh 


L  Z 

X  X 


+  L  Z 

y  y 


+  L  Z 
z  z 


For  the  residual  of  ^  A  h: 

^  ^  ~  (^ELl  A<|^  +  eL2  cosc^  ATv  +  eL3  A 
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or 


where 


and 


CE(1) 

= 

ELI 

CE(2) 

=  ^ 

EL2 

CE(3) 

= 

EL3 

RES 

II 

'Ah 

ELI 

=  CK3  (dX 
xh 

sin  <|) 

EL  2 

=  -(C  +  H) 

®yh 

EL3 

/-Sj 

=  -  °zh 

T 

^xh 

yh  zh 

u  = 
yh 

-  A  ,  L  . 
xh  zh 

11 

iJ 

1 

zh 


-  CK4  cos  (|)  sin(^) 


For  the  residual  of  A 


or 


where 


r 


RDl  A<|)  +  RD2  cos  ^  A>  +  RD3  A 


CE(1)  =  Vw^ 

RDl 

CE(2)  = 

RD2 

CE(3)  =  /W^ 

RD3 

RES  = 

t> 

RDl  =  -  CK3  cos©  ( 

RDl  =  -  CK3  cos©  ^ /^y  ®  +  x  -  jOL^,^  +  CK3  sin© 


+  CK4 
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RD2  =  (C  +  H)  Qos  ©  "  y  +  /5Ly )  +  sin  0  (  0  +  x  -  jOL^) 

RD3  =  cos  (|)  Q  cos  0  (^y  ©  ^  ®  ^  f'x®  ’  ^ 

+  sin  (j?  p-  z  +  p 

XI  =  -  Y1  0 
Y1  =  XI  © 

(^x  =  X  +  XI 
Py  =  y  +  Y1 

Pz-  ' 


c .  The  Least  Squares  Matrix 

The  least  squares  matrix,  Q,  has  the  form: 


Q(l) 

Q(2) 

Q(3) 

Q(2) 

Q(4) 

Q(5) 

Q(3) 

Q(5) 

Q(6) 

The  normal  equations  are: 


P(l) 

P(2) 

=  &] 

AcJ) 

cos^A^ 

P(3) 

Ah 
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These  equations  are  solved  for  1^^°®  4*  ^  H  by  pre¬ 

multiplying  both  sides  of  the  equation  by[Qj 

In  case  one,  two,  or  all  three  of  the  original  calculated  cor¬ 
rections,  /i<^and  coscj)  A>  ,  or  A  H  are  greater  than  a  specified  limit, 
the  logical  procedure  is  as  follows:  Take,  as  an  example,  the  case  when 
A  H  fails  to  pass  the  test.  That  is*  limA<^,  cos<j3  <  lim 

cos<{)A?v  ,  and  AH  >  lim  A  H.  In  this  case  we  replace  the  magnitude  of 
A  H  with  its  limit  value,  lim  A  H,  and  proceed  to  solve  for  the  remaining 
corrections . 


From  the  complete  least  squares  formulation  we  have: 


P(l) 

Q(l)  Q(2)  Q(3) 

A  4 

P(2) 

= 

Q(2)  Q(4)  Q(5) 

cos  a\ 

P(3) 

_Q(3)  Q(5)  Q(6)_ 

A  H 

Now  replace  the  quantity  A  H  by  its  limit  and  the  sign  of  AH 
and  then  rearrange: 


P(l) 

ah 

■  1^H| 

lim 

AH  Q(3)  =  Q(l)  A^) 

+  Q(2) 

cos  AT' 

(11) 

P(2) 

AH 
■  lAH) 

lim 

AH  Q(5)  =  Q(2)  A4> 

+  Q(4) 

cos  A  T- 

(12) 

P(3) 

AH 

‘lAHl 

lim 

AH  Q(6)  =  Q(3)  A4> 

+  Q(5) 

cos^  A^ 

(13) 

Note  that  the  quantities  on  the  left  hand  side  are  known  and  that  only  A^) 
and  cosC^  AT'  remain  as  unknowns. 


The  last  of  the  series  of  three  equations,  equation  (13),  can  be 

neglected  . 


This  can  be  verified  by  repeating  the  least  squares  procedure  for 
two  unknown  . 

There  are  now  two  equations  in  two  unknowns  --  (Jl)  and  (12)-- 
which  have  tbe  form 
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P(l) 

Q(l) 

Q(2) 

a4> 

P(2) 

Q(2) 

Q(4)_ 

cos(j) 

P(l) 

A(l) 

A(2) 

A 

cos(|) 

=  [*] 

a4> 

P(2) 

A(3) 

A(4)_ 

cos(|  A% 

where 

A(l)  =  Q(l) 

A(2)  =  Q(2) 

A(3)  =  Q(2) 

A(4)  =  Q(4) 

from  which  it  is  possible  to  solve  for  A  4^  and  coscpA^ 


Ac() 

cos  4? 


P(l) 

P(2) 


The  values  of  and  cos  <J>  A>  obtained  in  this  manner  are  the  best  values, 

in  a  least  squares  sense,  that  accommodate  the  limiting  value  of  AH.  This 
same  procedure  would  be  followed  should  either  A ^  or  cos  fail  the 

limit  check. 


What  happens  if  two  of  the  calculated  quantities  Ad}  ,  COS  dj)  A>. , 
or  A  H  fail  to  pass  the  limiting  test?  Again,  assume  AH  to  have  failed 
but  also  assume  A4>  to  have  failed. 


The  least  squares  equations  become: 

ACt  1  Ah 

P(1)  -  Q(l)  M  -  Q(3)  ah 

^(2)  -  Q(2)  jl^jlim  A4>  -  Q(5)j-|||lim  Ah 

||f  I  ^  ” 


Q(2)  cos^A7^(14) 
Q(4)  cosit)A>i(15) 
Q(5)  cos  4)A>(16) 
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These  equations  are  all  equivalent  to  the  least  squares  equation 
for  finding  one  unknown  from  a  group  of  measures,  equation  (15).  Finally 
we  have 


cos 


c|) 


-  o<^>  illl 


11m  Acjy  -  Q(5),^^llm  A  H 
Q(4) 


which  Is  a  weighted  average  for  cos  A  7^  ,  given  the  limiting  values  of' 
A4  and  AH.  A  similar  procedure  Is  followed  regardless  of  which  two  fall 
the  test. 

If  all  three  of  the  original  corrections  fall  the  test,  then: 

=  rcoi|"l>\ 

If  one  or  two  of  the  original  calculated  corrections  fail  to  pass 
the  test  and  new  least  square  values  are  calculated  for  the  quantities, 
these  new  values  are  also  checked  against  the  limiting  values  and,  if  they 
fail,  a  procedure  similar  to  that  just  described  is  followed.  The  flow 
diagram  explains  what  happens  if  the  newly  calculated  values  fail  the  test 
in  a  reasonable  amount  of  detail. 


d .  Corrections  to  the  Geocentric  Station  Location  Vector: 

The  corrections  to  the  geocentric  station  location  vector,  for 
station  L,  are  calculated  as  follows: 


.Ax  (L) 
s 

Ay^(i) 

Ae^(l) 


(C  +  H)  stn  >  (S2)  -  cos  X  [^1)  sin  (SI)  +  cos  C|)  (53^ 

-  (C  +  H)  cos  A  (S2)  -  sin  a[(CK1)  sin  (SI)  +  cos  (53^ 

-  (CK2)  cos  (p  (SI)  +  sin  (j)  (53) 
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where : 


SI  =  A<()  or 


a4> 


lim  A  4* 


52  =  cosc^  or  I 

A  „  ah  ,  . 

53  =  AH  or,— 

AH| 


cos(i>  AX 


cos  <p  A>\ 
AH 


lim  cosc^  A'?v 


The  values  of  AX  (L)  ,  AY  (L)  and  A  2  (L)  are  obviously  governed  by  the 
value  of  SI,  S2 ,  or  S3  selected.  The  proper  value  of  SI,  S2 ,  or  S3  is  in 
turn  governed  by  the  checking  procedure  used  on  the  original  calculated 
values  of  Ac^  ,  cosi?  A>  ,  and  AH.  This  procedure  is  illustrated  in  the 
flow  diagram  and  in  the  preceding  subsection  (c) . 


e .  Calculation  of  Final  Station  Coordinates 


The  final  values  of  the  station  coordinates  are  calculated  as 

follows; 

4^  =  Cfj  +  SI 

^  =  A  +  S2/cos  (!|> 

H  =  H  4  S3 

where  SI,  S2 ,  and  S3  are  the  same  as  those  described  in  (d)  above. 
The  manner  in  which  the  proper  value  of  SI,  S2,  or  S3  is  selected  is  identical 
to  the  method  explained  there. 


f .  Datum  Translation  Vector 


The  datum  translation  vector  is  calculated  in  the  following  manner: 

2,  ■  "x  <>•> 

■  "5  W  (1.) 

•^L  S 
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y  AYgCL)  •  Wy  (L) 

Ay  =  — - 

5  W  (L) 

L  S 


)>  AZg(L)  •  (L) 

L  S 

V  W  (L) 

^  L  S 


where  W  ,  W  ,  and  W  are  the  weighting  factors  for  each  component  of 
the  station  location  correction  vector  (i.e.,  AXg.  AYg,  AZg). 


Correction  to  Station  Location  Vector 


The  correction  to  each  individual  station,  to  determine  its  loca¬ 
tion  with  respect  to  the  datum,  is  calculated  as  follows: 


A  Xg/^  (L) 


=  AXg(L)  - 


AXp  (L) 


where  the  subscript  S/D  implies  "station  with  respect  to  the  datum"  and, 
as  before,  the  subscripts  S  and  D  imply  station  and  datum  where  both  are 
referenced  to  a  geocentric  coordinate  system. 


3.3  FLOW  DIAGRAMS 

The  following  charts  are  similar  to  those  of  Section  2.3.  Since 
the  Station  Locator  Program  was  written  in  FORTRAN  and  ALTAC .  the  symbolic 
program  location  names  used  in  the  connecters  are  numbers  rather  than  let¬ 
ters.  Arrows  at  the  bottom  of  one  page  Indicate  that  the  flow  is  continued 
as  indicated  by  the  corresponding  arrow  at  the  top  of  the  next  page. 
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or  THE 
QUANTITIES 


KLPLC  sO 

V  f  ^ 


WoJenO  of  -PROCKAn 


CALL  EXIT 


INITIALIZATION  Or 
LEAST  SQUARES  MATRIX 
CLEMENTS  4  LOGIC  PLOW 
PARAMETERS . 

■pto-o 

PCl)»0 

PCS)  s  o 

Q<0«O 

Qcg«o 

Q.o)-o 


QC4)«o 

Q<5)-0 

<a<;«)»o 
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Yes 
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STATION  ANO  DATUM  INPOAMATION 
DATUM  lO^  STATION  ID, 

Li'wt  ON  j  Li»*y  ON 

Lim  on  ah,  ^the  weighting 
factors  Wx«  ,  Wy, 4  W>, 


^NO 


NCK 


■■o 


IS 

h)  • 


NO 


IS 


ESTAVUSH  DATUM 

g  ^^CNSE 

check; 

Last*  DATUM  ID 

NO 
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NO 
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_ i _ 

CALCULATE  THE  QUAt-TTlTIES 
CE  o>,  ceco,  ceC3\  4  "RES  RsR 
THE  RESIDUAL.  />AK _ 


ENTER  SUBROUTINE  SUB  1 
RETURN  EROM  SUB  1 


CHECK  THE  WEI6HTIN6  TACTOK 
ON^.  XT  ZERO,  SET  fT  EQUAL 
TO  1  .CALCULATE  THE  QUANTITIES 
CE  0) ,  CCCZ>,  CE  Cl> ,  4  -RES  TOR 
THE  RCSIOUAL 


ENTER  JUBRPtfTINE  L 
RETURN  FROM  SUB  1 


calculate  the  value  of  the 
determinant  of  the  least 

daUARES  3x3  MATRIX  -  DCT 
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NO 


^YC5 
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IShT 
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[ShT 

AO)- <3(0  ,  A<z)-aa) 
a(3)-<3.w  ,  A  CO  -  gw 
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ao) 


Q.C**) 
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5(1) 

5C3)»  .  IJm.  AH 

lAHl 


^asoo^ 


S<2)  »  cos  ^  APt 
5 a)  «  i^^.lim  A^ 
5(3)  =,  AH  Li^  AH 

(StTr _ 


^8600 


SCI)  «  A9 

/  «o 

S  (Z)  m  COS^AA 

"\8700 

50)»  zlH 

CALCULATE  AK»(0^  AYsCO^  AZ»Cl). 
THESE  A  AC  THE  CORRECnONS  TO 
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3.4  INPUT  AND  OUTPUT  FORMATS 
Residual  Information 


Columns 


Quantity  and  Explanation 


Format 


1-4 


Fixed  Point  Constant  (FLAG) 

0  implies  o(  ,  ^  residuals 
+  N  Implies  A,  h  residuals 
-  N  implies  ja  residuals 


I  4 


5-8 

Station  Identification  Number 
(fixed  point  constant) 

I 

4 

9-21 

Time 

(min)  -  Time  of  Observation 

E 

13.8 

22-34 

AO 

• 

or  -  (Km  or  KM/sec) 

E 

13.8 

35-47 

^cos  h  A  A  or  pcos$Ao<-  (Km) 

E 

13.8 

48-60 

/Ah 

or  -  (Km) 

E 

13.8 

61-64 

Wp 

or  W  . 

P 

F 

4.2 

65-68 

“a 

or  M  ^ 

F 

4.2 

69-72 

or 

F 

4.2 

K(.:  (' 


Range  residuals  can  be  processed  with  either  type  of  angles  o'  ^ 
or  A,  h;  range-rate  information  must  be  processed  separately. 


Note  2 


Residuals  and  weighting  factors  must  be  entered  in  the  designated 
locations.  If  a  residual  does  not  exist  for  a  particular  case,  the  corres¬ 
ponding  location  should  be  blank. 
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General  Information  -  Input 


Columns 

(Quantity  and  Explanation 

Format 

1-14 

Theta  Greenwich  (6gr)  -  the 

Greenwich  sidereal  time  at  epoch 
time  tQ  (degrees) 

F  14.0* 

15-72 

Blank 

♦Recall 
format . 

that  a  punched  decimal  will  override  the  decimal  stated 

in  the  F 

Station 

Information  -  Input 

Columns 

Quantity  and  Explanation 

Format 

1-4 

Datum  identification  number 

I  4 

5-8 

Station  identification  number 

I  4 

9-17 

Latitude  of  station  (J  -  degrees 

F  9.0 

18-27 

Longitude  of  station  ^  g  ■  degrees 

F  10.0 

28-33 

Height  of  station  H  -  meters 

F  6.0 

34-40 

Limit  of  e  -  degrees 

F  7.0 

41-48 

Limit  of  cos  -  degrees 

F  8.0 

49-54 

Limit  of  ^H  -  meters 

F  6.0 

55-60 

Weighting  Factor  for  X  component  of 
station  location  vector  -  Wxg 

F  6.0 
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Station  Information  -  Input  (continued) 


Columns 


61-66 


67-72 


Quantity  and  Explanation 

Weighting  factor  for  Y  component  of 
station  location  vector  - 

s 

Weighting  factor  for  Z  component  f 
station  location  vector  -  W^ 

s 


Format 


F  6.0 


F  6.0 


Position  and  Velocity  Data 

The  ephemeris  data  is  input  to  this  program  as  a  binary  tape.  This 
tape  is  output  in  the  proper  format  from  the  Ephemeris  Program  of  Task  1.  Its 
format  is  as  follows:  time  of  observation,  position  of  satellite,  and  velocity 
of  satellite.  Position  and  velocity  are  given  in  a  geocentric  coordinate  sys¬ 
tem  with  units  of  Earth  radii  and  Earth  radii  per  kemin.  Time  is  given  in  min¬ 
utes.  In  symbolic  form,  the  input  format  is  t,  x,  y,  z,  x,  y,  z. 


The  position  and  velocity  data  is  input  on  a  binary  tape  through 
drive  number  3. 


The  quantities  --  residual  information,  general  information,  and 
station  information--  are  input  into  the  present  program  from  tape  unit  num¬ 
ber  0.  This  tape  can  be  prepared  from  punched  cards  of  the  formats  listed 
above.  Their  order  of  input  --  more  than  one  datum  and  station  --  is  as  fol¬ 
lows  : 
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3LANK 

BLANK 


(  -RES.  INFO. 
STATION  L 
{  OTA.  INFO. 
STA.  L  JATUM  N 

BLAN  K 


/  KES.  INFO. 
STATION  K 
/  STA.  INFO. 
STA.K  JATUM  N  - 

blank 


(  RES.  INFO. 

I  STATION  J 
{  STA.  INFO. 

I  STA.T  DATUM  M  - 

blank  1 


i 


{ - ftES.j;NFa 

STATION  I 
'  STA.  Info. 

STA.  I  DATUM  M  _ 
GEN.  INFO.  1 


The  stations  of  a  particular  datum  do  not  have  to  be  in  numerical 
order;  however,  all  the  observation  residuals  of  a  station  and  all  of  the 
various  stations  of  a  datum  must  be  sorted  and  input  immediately  after  the 
corresponding  station  and  datum  information. 
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The  quantities  in  the  output  are  as  follows; 

Page  Output  for  Each  Station 

1  Datum  ID  .  NDATID 

Station  ID  ....  NSTID 

Assumed  station  latitude  (degrees)  .  Cf) 

Assumed  station  longitude  (degrees)  .  "X 

Assumed  station  height  (meters)  .  H 

Corrected  station  latitude  (degrees)  ....  4^ 

Corrected  station  longitude  (degrees)  ... 

Corrected  station  height  (meters)  .  H 

Correction  to  assumed  station  latitude  (degrees)  . . 

Correction  to  assumed  station  longitude  (degrees)  . 

Correction  to  assumed  station  height  (meters)  . AH 

Limit  on  correction  to  assumed  station  latitude  (degrees)  . . .  lim  Adj 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 
17 


cos  (|;  /cos  cp 


Limit  on  correction  to  assumed  station  longitude  (degrees)  ..  lim  cosC|? 

/cos  ^ 

Limit  on  correction  to  assumed  station  height  (meters)  . lim  A  H 

Final  station  latitude  (degrees)  . cj? 

Final  station  longitude  (degrees)  . 

Final  station  height  (meters)  .  H 


Output  for  each  Datum 

13  The  assumed  rectangular  coordinates  of  each  station  and  their  corrections 
(megameters)  -  X  (L) ,  Y  (L) ,  Z  (L) ,  /iX  (L),  A Y  (L) ,  and  A  Z  (L) 


19 


Datum  correction  vector  (megameters)  ....  Ax^,  AY^,  ^ 
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Output  for  each  Datum  (continued) 


20  Station  translation  vector  with  respect  to  the  datum  for  each  station, 
(megameters)  ^  X  ,  AY^,  AZ  . 

S  S  5 

This  output  is  obtained  by  printing  tape  5  with  data  select  0. 

3.5  OPERATING  PROCEDURE 

The  deck  submitted  at  execution  time  has  the  following  make  up; 
Starting  in 


Col.  17 

Col.  25 

Purpose 

JOB 

xxxxx. . . 

Signals  the  start  of  a  new  job.  Cols. 25 
on  are  typed  on  the  flexowriter  at  the 
beginning  of  the  program  and  may  contain 
any  desired  alphanumeric  information. 

REWIND 

0,3 

Rewinds  tapes  0  and  3. 

RPL 

1,  DATA, GO 

Reads  routine  named  DATA  from  tape  1 
into  memory  and  transfers  control  to  it. 

TAPE 

0 

Control  instruction  for  DATA  telling  it 
which  tape  to  put  the  following  data  on. 

.... 

1 

1 

1 

J 

Data  to  be  put  on  tape  0. 

ENDDATA 

Control  instruction  for  DATA  telling  it 
that  there  is  no  more  data  to  put  on 
tape  0. 

REWIND 

0 

Rewinds  tape  0 

RPL 

4, XXX. . .XXX.GO 

XXX... XXX  is  the  identity  of  the  Station 
Locator  Program  (up  to  16  characters) .  It 
is  read  into  memory  from  tape  4  and  control 
is  transferred  to  it. 
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Tapes  Used; 


Logical  Tape  No. 


Data  input. 

Binary  input  tape  containing  t-t  ,  r,  and  r  at 
observation  times  (output  from  Orbit  Correction 
Program) . 

RPL  tape  containing  Station  Locator  Program. 
Output  tape. 


The  use  of  the  Station  Locator  Program,  together  with  the  Orbit 
Correction  Program  is  outlined  in  Section  2.5. 


3.6  EXPERIMENTATION 

Two  test  cases  were  used  to  check  the  use  of  the  Station  Locator 
Program,  in  conjunction  with  the  Orbit  Correction  Program. 

(1)  Using  the  Differential  Correction  Program,  observations  were 
simulated  from  two  stations  for  the  following  orbit: 


L  =4.890  210  9 
o 

a  „  =  -.007  371  516  8 


h  =  -  .618  213  08 

X 

o 

h  =  .484  673  16 


yN  =  .007  750  635  7  1 

^  o 

t  =  0.0  minutes,  98  days,  1961 


,704  418  46 


The  coordinates  of  the  two  stations  were: 


Station  44 

C|)  =  27?020  3 
^g=  279?886  9. 

=  12  meters 


Station  3 


4^  =  35?707  2 
=  139?491  7 


H  =0  meters 
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Using  these  observations  and  the  same  orbital  elements,  the  Orbit  Correc¬ 
tion  Program  was  used  to  compute  residuals  by  altering  the  station  coordi¬ 
nates.  The  residuals  generated  were  four  slant  ranges,  four  pairs  of  right 
ascension  and  declination  and  four  pairs  of  azimuth  and  elevation  angle  for 
Station  3  and  three  slant  ranges,  three  pairs  of  right  ascension  and  decli¬ 
nation  and  four  pairs  of  azimuth  and  elevation  angle  for  Station  44.  The 
root  mean  square  of  the  residuals  in  these  observations  caused  by  "moving" 
the  stations  was  0.1714  km.  The  amount  of  the  station  displacements  is 
shown  in  Table  1  below. 

These  residuals  were  used  in  the  Station  Locator  Program  to  cor¬ 
rect  the  station  coordinates.  The  limits  on  the  correction  were  set  at 
very  large  numbers  so  that  all  three  coordinates  would  be  corrected  for 
both  stations.  These  corrected  station  coordinates  were  fed  back  into  the 
Orbit  Correction  Program  to  get  new  residuals  which  were  again  fed  into 
the  Station  Locator  Program  to  further  correct  the  station  locations.  The 
results  are  summarized  in  the  following  table. 


TABLE  1 

RESULTS  OF  EXPERIMENTATION 


Displaced 

Station 

Station  44 


)  Acll| 

-  degrees 

-  meters 

27.1 

166 

D21  8 
.7 

^  E 

-  degrees 

279 

.885  4 

1 

-  meters 

148 

.5 

H 

-  meters 

112 

.0 

1  A  Hi 

-  meters 

100 

.0 

Stat ion 

3 

-  degrees 

35. 

705  2 

-  meters 

222 

.2 

>E 

-  degrees 

139 

.493  2 

1  A>1 

-  meters 

135 

.4 

H 

-  meters 

200 

.0 

1  Ah  1 

-  meters 

200 

.0 

After  First 

After  Second 

Correction 

Correction 

27.020  241 

27.020  299 

6.6 

.1 

279.887  01 

279.886  90 

10.9 

0 

17.109  621 

12.266 

5.1 

.3 

35.706 

880 

35.707  172 

35.6 

3.1 

139.491 

62 

139.491  69 

7.2 

.9 

27.791 

097 

-  6.695 

27.8 

6.7 
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(2)  Using  the  same  simulated  observations  as  before  and  the  same 
displaced  station  coordinates  as  before,  but  with  the  following  orbital 
elements; 


=  4.891 

h 

=  -  .618 

2 

0 

X 

xN 

=  -.007 

37 

h  ° 
y 

=  .484 

7 

o 

0 

yN 
^  o 

=  .007 

75 

h 

=  .704 

4 

z 

0 

four  differential  corrections  were  made  with  the  Orbit  Correction  Progra*. 
to  correct  the  orbital  elements  as  much  as  possible  before  correcting  the 
station  locations.  The  final  orbital  elements  were: 


0 

=  4.890 

224 

84 

h 

X 

0 

00 

1 

183 

83 

xN 

o 

=  -.007 

377 

960 

h 

It 

00 

681 

51 

yN 

o 

=  .007 

746 

031 

h 

z 

0 

■=  .704 

00 

41 

The  initial  root  mean  square  of  the  range  and  angle  residuals  was  20.377  km 
and  for  the  range-rate  residuals  it  was  .030  982  km/sec,  while  the  final 
root  mean  square  for  the  range  and  angle  residuals  was  .141  km  and  for  the 
range-rate  residuals  it  was  .000  177  km/sec.  The  final  residuals  were  fed 
into  the  Station  Locator  Program  to  correct  the  station  coordinates.  Four 
observations  of  range  rate  from  each  station  were  used,  in  addition  to  the 
same  observations  used  in  the  above  experimentation.  The  results  are  sum¬ 
marized  in  Table  2. 
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TABLE  2 

RESULTS  OF  EXPERIMENTATION 


Before  Station  After  Station 

Correction  Correction 


Station 

44 

1  A45\ 

-  degrees 

27.021  8 

27.020  296 

-  meters 

166.7 

.4 

E 

-  degrees 

279.885  4 

279.887  01 

1  A'M 

-  meters 

148.5 

10.9 

H 

-  meters 

112.0 

19.247  218 

1  Ah  I 

-  meters 

100.0 

7.2 

Station 

3 

. 

-  degrees 

35.705  2 

35.706  857 

-  meters 

222.2 

38.1 

-  degrees 

139.493  2 

139.491  71 

(  A9<l 

-  meters 

135.4 

.9 

H 

-  meters 

200.0 

-  32.768  020 

IAh  I 

-  meters 

200.0 

32.8 
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SECTION  4 

MODIFICATIONS 


This  section  is  devoted  to  suggested  modifications  to  the  programs 
described  above  and  to  subsidiary  progrfims  which  may  prove  desirable. 

The  programs  heretofore  described  use  a  single  Earth  ellipsoid  for 
all  datums.  The  station  coordinates  produced  by  the  Station  Locator  Program 
will  be  correct  in  cartesian  coordinates  and  consistent  with  that  ellipsoid 
in  latitude,  longitude  and  height.  All  station  coordinates  should,  in  the 
near  future,  be  recomputed  on  the  basis  of  such  a  world-wide  datum  or,  alter¬ 
native,  available  in  a  cartesian  system. 

There  exist,  however,  good  observations  (e.g.,  surveys)  which  are 
not  used  by  the  present  programs.  They  connect  many  more  stations  than  can 
observe  the  geodetic  satellite.  The  presence  of  such  ties  is  felt  in  the 
Station  Locator  Program  in  the  limits  imposed  on  the  station  coordinate  cor¬ 
rections.  An  alternative  method  of  accounting  for  such  additional  informa¬ 
tion  as  the  geodetic  ties  is  to  represent  it  by  equations  of  condition  which 
state  that  each  station  is  where  the  surveys  put  it.  Initially  these  equa¬ 
tions  will  say  for  every  station  coordinate; 

AX.  =  0  . 

J 

Each  equation  should  be  weighted  according  to  the  accuracy  of  the  previous 
determination  of  that  coordinate.  After  the  first  correction,  the  right 
sides  of  the  equations  would  no  longer  be  zero.  Each  equation  then  would 
be  equal  to  the  negative  sum  of  all  the  previous  corrections  to  that  coor¬ 
dinate.  The  advantages  of  this  procedure  include: 
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(1)  No  station  position  is  ever  underdetermined. 

(2)  Because  the  limits  are  not  used,  there  is  but 
one  solution. 

(3)  A  new  uncertainty  can  be  generated  on  the  basis 
of  both  the  old  ties  and  the  new  geodetic  satel¬ 
lite  observations. 

The  third  feature  can  be  added  to  the  present  Station  Locator 
Program  only  insofar  as  the  new  observations  are  concerned.  If  valid 
estimates  of  the  uncertainties  in  each  observation  are  used,  the  uncertainty 
in  each  station  coordinate  can  be  generated. 

The  uncertainties  in  the  observations  can  also  be  translated  into 
uncertainties  in  orbital  elements.  This  feature  can  be  added  to  the  Orbit 
Correction  Program.  At  the  same  time  the  observational  uncertainties  can 
be  used  to  weight  the  corresponding  equations  of  the  Orbit  Correction  Pro¬ 
gram.  These  feat  ires  will  improve  the  solution  for  the  elements,  give  each 
element  a  figure  of  merit  and  enable  the  program  to  be  used  to  evaluate 
various  geodetic  satellite  orbiiswith  respect  to  their  accuracy  in  locating 
a  given  station. 

If  these  weighting  factors  are  introduced  into  the  Orbit  Correc¬ 
tion  Program, the  criteria  for  the  rejection  of  discordant  observations  should 
be  re-examined.  Rejection  and/or  weighting  on  the  baSiS  of  the  distribution 
and  numbers  of  observations  should  be  considered. 

Other  changes  could  be  introduced  into  the  satellite  orbit  model 
to  improve  accuracy.  These  include  the  introduction  of  the  perturbative 
acceleration  due  to  high-order  zonal  and  sectorial  harmonics  of  the  Earth's 
gravitational  potential  and  due  to  the  pressure  of  solar  radiation,  as  well 
as  refinements  in  the  computation  of  aerodynamic  drag.  The  effects  of  lunar 
and  solar  attraction  and  of  magnetic  forces  should  I)e  evaluated  and  intro¬ 
duced  if  necessary. 

As  indicated  in  the  Introduction,  the  two  programs  should  be  com¬ 
bined  in  a  flexible  manner  so  that  only  orbital  elements  or  orbital  elements 
and  station  coordinates  simultaneously  could  be  corrected. 

The  values  of  the  coefficients  of  the  harmonics  of  the  Earth's  gravi¬ 
tational  potential  may  be  determined  by  differential  correction  processes  also. 
Whether  this  should  be  done  simultaneously  with  other  corrections  is  subject 
to  the  same  sort  of  argument  as  presented  in  the  Introduction. 
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Appendix  A 


Glossary  of  Symbols,  Constants,  and  Units 


Definition 


Mean  distance  or  semi-major  axis  of 
elliptical  orbit 

Earth's  equatorial  radius  =  6378.150 
meters 

Vector  directed  to  perigee  having  magni¬ 
tude  e,  and  with  orbit  plane  components: 
axn*  along  N;  and  ay^.  along  M 

Azimuth;  measured  in  a  positive  sense  from 
the  north  point  eastward  in  the  horizon  plane 

Unit  vector  in  the  equator  system  forming 
an  orthogonal  set  with  L  and  D;  components 
are  Axi  Ay,  A^  in  the  K  coordinate 

system 

Unit  vector  in  the  horizon  system  forming 
an  orthogonal  set  with  L  and  tSh;  components 
are  Axj^.  A^j^  in  the  S,  ^  coordinate 

system 

Atmospheric  drag  coefficient 

Reference  atmospheric  drag  coefficient  = 
0.92 

Diameter  of  satellite 

Unit  vector  in  the  equator  system  forming 
an  orthogonal  set  with  L  and  A;  components 
are  D^.  Dy,  in  the^,  ^  system 

Unit  vector  in  the  horizon  system  forming 
an  orthogonal  set  with  L  and  components 
are  Dyj^.  in  the  ^  system 

Eccentricity  of  elliptical  orbit;  also  used  as 
the  earth's  meridional  eccentricity, 
e^  =  2f  -  f^  where  f  is  the  flattening 


Units 

earth  radii 

earth  radii 

dimensionless 

radians 

dimensionless 

dimensionless 

dimensionless 

dimensionless 

meters 

dimensionless 

dimensionless 

dimensionless 


Symbol 


E 

E 

f 

h 

h 

H 


H' 


i 

I 


J' 


J 


5 


J 


k 


e 


K' 


Definition 


Eccentric  anomaly 

Unit  vector  in  horizon  system  directed 
toward  east 

Flattening  of  the  spheriod  (Earth)  =  1 

19873 

Elevation  angle;  angular  distance  above 
(positive)  or  below  (negative)  the  horizon 

A  vector  related  to  the  angular  momentum 
(per  unit  mass),  being  equal  to  (£  x  £)/j^5 — 

Height  above  sea  level,  or  height  above  the 
geoid 

A  coefficient  of  the  3rd  order,  describing 
the  perturbative  acceleration  resulting  from 
the  3rd  harmonic  of  the  Earth's  potential  = 
-.00000575 

Inclination  of  the  orbit  plane  to  the  plane  of 
the  equator 

Unit  vector  in  the  inertial  framework 
directed  toward  the  vernal  equinox 

Coefficient  of  the  2nd  order  zonal  harmonic 
of  the  Earth's  potential  =  .00162341 

Coefficient  of  the  5th  order  harmonic  of  the 
Earth's  potential  =  -.0000002 

Unit  vector  in  the  inertial  framework  in  the 
equator  plane  forming  an  orthogonal  set  with 
and  ^ 

Geocentric  gravitational  constant  = 

.07436574 

Coefficient  of  the  4th  order  harmonic 
involving  "flattening"  of  spheroid  and  centrif¬ 
ugal  acceleration  =  .00000795 


Units 


radians 

dimensionless 


dimensionless 


radians 

radii  ^/k’^inm^^^ 


earth  radii 


dimensionless 


radians 


dimensionless 


dimensionless 


dimensionless 


dimensionless 


...  3/2.,,  -1  .  , 
radii  /(kg  min) 


dimensionless 
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Symbol 


Definition 


Units 


Unit  vector  in  the  inertial  freumework 
directed  northward  along  the  polar  axis 
of  the  Earth 

dimensionless 

True  orbital  longitude,  referred  to  intertial 
coordinate  system  1  = 

radians 

Mean  orbital  longitude  referred  to  inertial 
coordinate  system  L  =  M  +  + 

radians 

Unit  vector  directed  from  observer  to 
to  satellite 

dimensionless 

Mass  of  orbiting  vehicle 

kilograms 

Mean  anomaly 

radians 

Molecular  weight  of  atmosphere 

Unit  vector  in  orbit  plane  directed  to  a 
point  90°  from  node  in  direction  of  motion; 
components  M^,  M  ,  in  the_I,  J,  K 
system 


gram  mol.  wt. 


dimensionless 


Mean  angular  motion 


radians  /min 


Unit  vector  to  ascending  node;  components  dimensionless 
Nx«  Ny,  in  the_I,  J,  K  system 

Semi-latus  rectum  earth  radii 

Unit  vector  directed  to  perigee;  components  dimensionless 
Px.  Py.  in  the^,  J^,  system 

Unit  vector  parallel  to  minor  Eixis  and  dimensionless 

velocity  vector  at  perigee;  components 
Qx>  Qy»  Qz  in  the_I,  £,  K  system 

Distance  of  object  from  geocenter  earth  radii 

Radial  component  of  vehicle  velocity  vector  earth  radii/k'^min 

Vector  from  geocenter  to  object;  components  earth  radii 
X,  y,  z  in  the  I,  J,  K  system  and  x^,  yw  in  the 
P,  Q  system 

Velocity  vector  of  object  relative  to  geocenter;  earth  radii/k  ^min 
components  x,  y,  z  in  the  I»  J,  K  system 
and  x^,  ^  in  the  P,  Q  system 
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Symbol 

Definition 

Units 

r'^ 

Total  perturbative  acceleration  of  vehicle; 
components  in  the  I,  J,  K  system 

earth  radii/(k*^ 
e 

min)' 

.\ 

Lb 

Perturbative  acceleration  of  vehicle  due  to 
Earth's  bulge;  components  Xg,  yg,  Zg  in 
the  I,  J,  K  system 

earth  radii/(k 

e 

min) 

.\ 

^JD 

Perturbative  acceleration  of  vehicle  due  to 
atmospheric  drag;  components  x^,  yp,  zq 
in  the  I,  J,  K  system 

earth  radii/(k^^ 

min) 

rh 

Normal  component  of  vehicle  velocity  vector 

earth  radii/k  ^min 
e 

rv 

Transverse  component  of  vehicle  velocity 
vector 

earth  radii/k  min 
e 

R 

Position  vector  of  geocenter  with  respect 
to  observer;  components  X,  Y,  Z  in  the  I, 

J,  K  system 

earth  radii 

k 

Velocity  vector  of  geocenter  with  respect 
to  observer;  components  X,  Y,  Z  in  the  I, 

J,  K  system 

,  -1 

earth  radii/kg  min 

s 

• 

Magnitude  of  r 

earth  radii/k  min 
e 

S 

Unit  vector  in  horizon  system  directed 

toward  south;  components  S„,  S  ,  S  in  the 

X  y  z 

I,  J,  K  system 

dimensionless 

t 

Time 

minutes 

T 
^  s 

Skin  temperature  of  satellite 

degrees  Kelvin 

u 

Argument  of  latitude;  u  =  v  +  A/ 

radians 

U 

Mean  argument  of  latitude;  U  =  M  + 

radians 

U 

Unit  vector  directed  along  radius  vector; 
components  Uy,  U^  in  the  I,  J,  K  system 

dimensionless 

V 

T  rue  anomaly 

radians 

^co 

Circular  satellite  velocity  at  unit  distance  = 

7. 9048  km/sec 

km/sec 
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Symbol 

Definition 

Units 

V 

Transverse  unit  vector  perpendicular  to 

dimensionless 

radius  vector;  components  V^,  V  ,  V  in 
the  I,  J,  K  system  ^  * 

W 

Unit  vector  perpendicular  to  orbit  plane; 

dimensionless 

components  V/  ,  W  ,  W  in  the  I,  J,  K 
system 

Z 

Unit  vector  in  horizon  system  directed 

dimensionless 

toward  zenith;  components  Z^,  Z^,  in 

the  I,  J,  K  system 

JU 

Right  ascension 

radians 

s 

Declination 

radians 

Emissivity  of  satellite  =  .9 

dimensionless 

& 

Local  sidereal  time 

radians 

Bgr 

Greenwich  sidereal  time 

radians 

& 

Rotational  rate  of  the  Earth  ^  .0043752689 

radians /min; 

A 

Longitude 

radians 

Mass  function  =  Vm^  +  m2 
(=  unity  for  m^^^m^) 

earth  mass 

Velocity  vector  relative  to  the  Earth's 
atmosphere;  components  Vy.,  m  the 

radii /k  min 
e 

I,  J,  K  system 

r 

Slant  range  to  vehicle 

earth  radii 

r. 

Atmospheric  density 

gm/cm 

t 

Range  rate  of  vehicle 

radii/k  ^min 
e 

(S 

Range  vector  to  vehicle;  components 

earth  radii 

I 

Px»  f^Y*  Pz,  ^  system 

• 

1  1 

-1 

p 

Velocity  vector  of  vehicle  relative  to 

radii/k  min 

r 

observer;  components  ix*  ^v’  P 

I,  J,  K  system  [  [ 

e 
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Symbol 

Definition 

Units 

O' 

Atmospheric  density  ratio  =  \  /fo 

where  ^  =  0.001225  gm/cm^  ' 

dimensionless 

cr 

s 

Stefan -boltzmann  constant  ■  5.672  x  10 

erg  cm'^deg"' 

Geodetic  latitude 

radians 

Jl 

Longitude  of  ascending  node  of  the  orbit; 
angle  measured  in  the  equator  plane  from 
the  vernal  equinox 

radians 

Argument  of  perigee 

radians 

The  above  symbols  are  modified  in  many  cases  by  the  following  symbols: 


Symbol 

Description 

Definition 

• 

Over  dot 

Denotes  time  derivative 
exclusive  of  perturbations 

\ 

Grave 

Denotes  time  derivative  due 
to  perturbations 

— 

Underscore 

Denotes  vector  quantity 

o 

Zero  subscript 

Denotes  quantity  referenced  to 
some  epoch  or  standard  value 

L 

Delta  prefix 

Denotes  residual  quantity 

c 

Subscript 

Denotes  computed  quantity 

obs 

Subscript 

Denotes  observed  quantity 
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Appendix  B 


List  of  Constants 


Orbit  Correction  Program 

and  Station 

.  Locator  Program 

Symbol 

Value 

Definition 

THGRO 

98. 67400833 

(I960)  ^ 

99.42093750 

(1961) 

►  ©GR  •  Greenwich  Sidereal  Time 
at  beginning  of  epoch  year 

99. 18221667 
98.94350000 

(1962) 

(1963) 

J 

SIDRT 

.9856472 

o  * 

0-360  ,  where  0  is  the  rotation  rate  of  the 
earth  in  degrees/mean  solar  day 

FLP  25 

-.25068448 

The  negative  of  the  rotation  rate  of  the  earth 
in  degrees/minute 

THDOT 

.05883447 

ThjB  rotation  rate  of  the  earth  in  radians/ 
kg 'min 

EMIS 

0.9 

The  emissivity  of  the  satellite  (used  in  the 
drag  calculation) 

XK 

.2504742  ElO 

A  constant  relating  the  units  used  in  the  drag 
calculation 

CDO 

.92 

An  emperically  determined  constant  used  to 
evaluate  the  drag  coefficient 

ONEPI 

1. 173913 

An  empirically  determined  constant  used  to 
evaluate  the  drag  coefficient 

SIXP9 

6.972  E9 

An  empirically  determined  constant  used  to 
evaulale  the  drag  coefficient 

CD 

2.0 

Initial  approximation  to  the  drag  coefficient 

SICS 

.5672  E-4 

the  Stefan-Boltzmann  constant 

XKE 

.07436574 

kg,  the  gravitational  constant 

F 

.0033523299 

f,  the  flattening  of  the  earth 

EPSQD 

6.6934216  E-3 

7  7 

e  (=2f-f  ),  where  e  is  the  eccentricity  of  the 
terrestrial  ellipsoid 
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Symbol 

Value 

Definition 

XJMPRM 

1.62341  E-3 

J',  the  2nd  harmonic  of  the  earth's  potential 

HMPRM 

-5.75  E-6 

H',  the  3rd  harmonic  of  the  earth's  potential 

XKMPRM 

7.96  E-6 

K',  the  4th  harmonic  of  the  earth's  potential 

XJAY5 

-0.2  e-6 

Jg,  the  5th  harmonic  of  the  earth's  potential 

XMPER'I 
XK2ER  ) 

.6378150  E7 

Meters  per  earth  radii 

XKS2RK 

7.9048 

km/sec  per  earth's  radii/ke  min 

VC03 

.49393823  E18 

3 

(V  )  where  V^^  is  the  speed  of  a  circular 
safe^llite  at  1  earth  radius,  in  (cm/sec) 

ERPM 

.  15678527  E-6 

Earth  radii  per  meter 

ER2MM 

.  15678527 

Earth  radii  per  megameter 

RHOO 

.001225 

P  ,  the  atmospheric  density  at  the  surface  of 

the  earth  in  gms/cm 


Reference  for  all  geodetic  constants:  Makemson,  Baker,  emd  Westrom  (AAS  Journal, 
Vol.  8,  #1,  p.l,  Spring  1961) 
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Appendix  C 


Denflity  Table 

(1959  ARDC  Model  Atmos. ) 

Used  in  the  Calculation  of  the  Drag  Perturbation  in  the 
Orbit  Correction  Prograjn 


H  (meters) 

Ln(r 

Me 

50.0  X  10^ 

-7.0310587 

28.966 

55.0 

-7. 6032281 

28.966 

60.0 

-8. 1538223 

28.966 

65.0 

-8. 7536457 

28.966 

70.  0 

-9.4124843 

28.966 

75.0 

-10. 143371 

28.966 

80.0 

-10.964450 

28.  97 

85.0 

-11.969229 

28.97 

90.  0 

-12.972557 

28.97 

95.0 

-14.027182 

28.94 

100.  0 

-15.003575 

28.90 

105.  0 

-15.911001 

28.  86 

110.0 

-16. 843617 

28.82 

115.  0 

-17.  625944 

28.77 

120.0 

-18.231593 

28.71 

125.  0 

-18.725451 

28.66 

130.0 

-19. 142311 

28.59 

135.0 

-19.503103 

28.53 

140.  0 

-19.820952 

28.45 

145.0 

-20. 125377 

28.  36 

150.  0 

-20. 361461 

28.27 

155.0 

-20.595472 

28.  16 

160.  0 

-20.802337 

28.04 

165.0 

-21.003537 

27.91 

170.  0 

-21.157902 

27.75 

175.0 

-21.305866 

27.57 

180.0 

-21.434534 

27.  36 

185.0 

-21.561085 

27.  12 

190.  0 

-21.685493 

26.  85 

195.0 

-21.807622 

26.59 

200.0  - 

-21.927794 

26.  32 

210. 0  X  10 

-22. 159273 

25.80 

220.0 

-22. 379078 

25.29 

235.  0 

-22.699344 

24.56 

250.0 

-23.008858 

23.87 

270.0 

-23.406146 

23.03 
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H  (meters) 

Ln  (/^ /A) 

Me. 

300.0 

-23.971776 

21.95 

330.0 

-24.504617 

21.06 

360.0 

-25.007847 

20.  33 

400.0 

-25.637420 

19.  56 

450.0 

-26.366557 

18.83 

500.0 

-27.040063 

18.28 

600.0 

-28.246917 

17.52 

900.0 

-31.867479 

16.50 

(/J>  =  0.001225  gms/cm^) 
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Appendix  D 


Error 


On-Line  Comment 


ERROR  IN  NODE 
COMPUTATION  0/0 

30  TIMES  THRU  LOOP 
WITHOUT  CLOSING  ON  E 


ERROR  IN  RUNGE  KUTTA 
ROUTINE 

PROGRAM  SHOULD  BE  ■ 
DUMPED 


END  OF  STATION  DATA 
CARDS 

ERROR  IN  TAN  A  DIVIDING 
0  BY  0 


STATION  NUMBER  XXXX 
NOT  FOUND 


NO.  OF  UNKNOWNS 
EXCEEDS  NO.  OF  OBSER¬ 
VATIONS 


Exits  and  On-Line  Comments* 


Origination 


Entered  from  XYZSB  if 
any  XNX  =  XNY  =  0 

Entered  from  XYZSB  if 
the  iteration  on  E+tU  does 
not  converge  to  within  10”^ 
after  30  iterations 

Entered  from  RUNGE  sub¬ 
routine  in  either  of  2  ways: 

a)  if  4t  *  0,  the  sign  of  the 

\  register  is  set  positive. 

b)  if  there  is  an  error  in 
computing  the  error  con¬ 
trol  term  (in  the  variable 
interval  mode),  the  sign 
of  the  A  register  is  set 
negative. 

Normal  exit  from  program 
if  simulation  is  being  used 

Entered  from  OBSIM  if: 

1.  XLX  =  XLY  =  0  (error 
in  computing  A),  or 

2.  XLSUBX  =  XLSUBY  =  0 
(error  in  computing  oi- ) 

Entered  from  ORBCO  if  the 
station  number  on  an  obser¬ 
vation  card  cannot  be  associ¬ 
ated  with  the  current  list  of 
station  numbers  in  the  pro¬ 
gram 

Entered  from  FUNCT  if  there 
are  less  good  (i.  e. ,  not 
rejected)  observations  than 
unknowns 


Action 


Normal  system  exit 
(location  4) 

Continues  with  30th 
value  of  E+tD 


Program  halts.  If 
advance  bar  is  pushed, 
program  transfers  to 
FINIS 


Program  transfers  to 
FINIS 

Program  transfers  to 
FINIS 


Reads  in  the  next  obser¬ 
vation  card 


Program  transfers  to 
FINIS 
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On-Line  Comment 


Origination 


Action 


(FINIS) 

ALL'S  WELL  THAT 
ENDS 


Entered  from  above 
referred  actions 


Clears  out  PUNCH  and 
PANT  buffers  and  executes 
the  normal  system  exit 
(location  4) 


*'^ius  table  applies  only  to  the  Orbit  Correction  Program;  the  Station  Locator 
Program  has  no  error  exits  or  on-line  comments. 
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